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Introduction 


Si. 


This paper is expanded from a mimeogra- 
phed memorandum circulated, and presented 
for discussion, at the spring meeting of the 
Seismological Society of Japan, at Tokyo, 
May 12, 1960. The intention is not to offer 
any new material, but to bring together more 
or less well known results in a form suitable 
for critical comment. 

During the remainder of his stay in Japan, 
the writer continued his study, reading, and 
discussion along the indicated lines. A few 
only of his personal indebtednesses are named 
in the section on Acknowledgments. 

The writer’s knowledge extends to only a 
small part of the vast body of seismological, 
geophysical, and geological data available in 
Japan. In spite of the great assistance he 
has received, he is certain that the paper 
“still contains inexactitudes, and even down- 
right mistakes. It is to be hoped that these 
are minor. Experienced workers will not be 
misled; but students are cautioned not to 
accept statements which they are inclined to 
~ question, without verifying them from the 
numerous publications which are more easily 
accessible to Japanese readers than to the pre- 
sent writer. 

The point of view is that of the writer’s 
text-book ‘‘Elementary Seismology’’; the pre- 
sent paper is in effect a revision of the cor- 
respoding discussion in Chapter 30 of the 
book. 

The presentation is with reference to the 
data of seismicity: accordingly, description 
and comparison refer to tectonic processes 
now going on, which in many areas differ 
significantly from those of the geological 
past. Moreover, present occurrence of earth- 


quakes in a given small area may be due to 
special circumstances (either geographically 
localized, or temporary in terms of the geo- 
logical time scale). Such local areas may 
present exceptions to the general tectonic pat- 
tern. 


§2. Nature of block and arc tectonics 


Most of the seismicity of the world is re- 
lated to two principal environments-those of 
block and arc tectonics; the latter is found 
chiefly in the circum-Pacific and Alpide belts. 
The work of Ewine and HEsgZzEN suggests that 
the seismicity of the Mid-Atlantic Ridge, and 
of other analogous oceanic ridges, may repre- 
sent a third principal tectonic type. 

From the viewpoint of historical geology, 
the distinction between regions of block 
tectonics and the active arc structures is usual- 
ly not one of essential character, but of the 
present stage of development. The materials 
of most active block structures have been 
subjected to successive orogenies; at one or 
more earlier periods such areas presumably 
presented all the characteristic features now 
associated with active arcs, suchas: (1) high 
seismicity including shallow, intermediate, 
and (in typical Pacific arc) deep earthquakes; 
(2) active volcanoes; (3) large gravity anoma- 
lies in narrow belts; (4) a pronounced foredeep. 
Block tectonics may be considered as a later 
stage which normally follows the folding and 
thrusting of a typical orogeny. This stage 
represents a lower degree of activity (especi- 
ally when conditions at the greater depths 
are considered) than the arc stage. Although 
large shallow earthquakes occur, and some of 
them break the surface in faulting, inter- 
mediate earthquakes are rare and deep earth- 
quakes absent (except when deep earthquakes 
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belonging to an adjacent active arc originate 
under the shallow block structure). Volcanoes 
are usually in a late stage of activity or ex- 
tinct, while foredeeps and gravity anomalies 
are less pronounced than in active arcs. Such 
criteria have been used in the following dis- 
cussion to distinguish the block tectonics of 
West Japan from the arc tectonics of adjacent 
areas. The distinction between block and arc 
tectonics is not perfectly exclusive; some 
seismic areas appear to be in a transitional 
condition, and others show characteristics of 
both kinds due to intersection or to complex 
geometrical relationship in three dimensions. 

Stratigraphy and other geological evidence 
in some instances suggests, and in others 
clearly demonstrates, that a present active 
arc formerly passed through a stage of block 
tectonics. 

This discussion has a definite relation to the 
concepts employed by the Soviet school of 
geologists and geophysicists. Active arcs are 
generally considered as representing active 
geosynclines. The major geosynclinal. belts 
are marginal to the stable platform areas. 
The regions of active block tectonics, which 
occur in the geosynclinal belts in juxtaposition 
to the active arcs, are classifiable as former 
platform areas, now involved in the develop- 
ment of geosynclinal movements, and being 
broken up in consequence. If this process 
goes far enough, the block structure may be 
completely destroyed, and the material recon- 
structed into a typical active arc. 


§3. Geographical subdivisions of Japan 


The exceptional importance of Japan in 
comparative tectonics derives from the pre- 
sence, in an accessible, civilized, and well- 
investigated region, of both block and arc 
structures, in such a geometrical and mechanic: 
al relationship as promises to shed light on 
the nature of both. 

There are five evident subdivisions: 

(1) West Japan, including the main islands 
west and southwest from the Fossa Magna; 
this is the major area of block tectonics, 
with only traces of arc features. 

(2) Northeast Japan, including Honshu east 
of the Fossa Magna, and southwestern Hok- 


kaido. This is part of a major Pacific arc. 
As in some other such arcs, there is block 
faulting in the interior; here it is chiefly in 
a belt adjacent to the Japan Sea. 

(3) Northeastern Hokkaido, which is the 
terminus of the active arc of the Kurile 
(Chishima) Islands. 

(4) The arc of the Ryukyu Islands, which 
extends into Kyushu, where it intersects and 
modifies the block structure. KOBAYASHI 
relates the structures of Korea centrally to 
this arc. 

(5) The Shichito arc, extending southward 
to the Ogasawara (Bonin) Islands. 

These divisions will now be discussed in- 
dividually, pointing out some of the more in- 
teresting comparisons with other regions. 
Later sections will sketch the block and arc 
tectonics of several other distant areas, pre- 
senting tentative general correlations with the 
Japanese region. 

In addition to the surface features of Japan, 
geophysical discussion must take into account 
the belts of deep and intermediate earthquakes 
discovered by Wapatl, as well as the structure 
of the crust as derived from seismological 
and gravity data. 

The study of most of the active Pacific 
arcs in hampered by their largely submarine 
character; often the accessible land area con- 
sists of a chain of small islands. This is true 
of the Shichito, Kurile, and Ryukyu arcs, 
except at their ends, where they intersect the 
other block and arc structures. Such areas 
of intersection are of special importance for 
investigation. There are three of them in 
Japan: one including the Fossa Magna and 
the western Kanto district, one in Kyushu, 
and one in Hokkaido. All three of these, 
Professor Tsusor points out, are associated 
with large negative gravity anomalies. The 
next similar critical areas in opposite direc- 
tions are in Kamchatka and Formosa (Tai- 
wan); findings there should be compared with 
Japan so far as possible. Kosayasut describes 
the transition from the Ryukyu arc to For- 
mosa rather as a sharp bend in the younger 
trends than as an intersection; whatever the 
interpretation, the contrast between adjacent 
block and arc structures is clear, 
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§4. West Japan 


In comparing West Japan with other areas 
of active block structure, a further subdivision 
must be considered. The longitudinal separa- 
tion, along the Median Tectonic Line, into an 
inner and an outer zone, is of very great im- 
portance in studying the geological past of 
West Japan; but it is not so significant in 
relation to present tectonic processes. On the 
other hand, there are important complications 
at the two ends of West Japan which do not 
appear in its central segment. At the west, 
although the block structures and the Median 
Line cross Kyushu, there is an evident in- 
tersection in three dimensions, where these 
shallow structures cross the Ryukyu arc. 
Among the known typical arc features are 
the clustering of shallow earthquakes, associat- 
ed with negative gravity anomalies, along 
the Bungo channel, and the central line of 
active volcanoes followed by epicenters of 
intermediate earthquakes. Continued study 
of the deep crustal structure of Kyushu, 
combining data from earthquake records, 
explosion seismology, and gravity observations 
on land and off shore, should clarify the geo- 
metrical relation of the intersecting structures. 
It is highly important to learn how the Mo- 
horovicic discontinuity behaves in such a 
region. 

A somewhat similar three-dimensional re- 
lation exists in the North Island of New 
Zealand, but the trends of block and arc 
structures are more nearly parallel there than 
~ in Kyushu. In the North Island, as in Kyu- 
shu, there is a central trough, largely filled 
with volcanic material; the active volcanoes 
and the epicenters of intermediate earthquakes 
follow the trough in both cases. One hypo- 
center is known at a depth of over 500 kilo- 
meters below the central part of the North 
Island; under Kyushu a few earthquakes are 
known at depths between 300 and 400 kilo- 
meters. 

In spite of some significant differences, the 
middle section of West Japan compares rather 
closely with other Pacific areas of dominant 
block tectonics, such as California and central 
New Zealand. The general level of small 


seismicity is not high; however, there are 
occasional very large earthquakes, especially 
off the southern coast. These earthquakes 
originate inshore from the Nankai Trench, 
which thus is comparable in position with the 
foredeeps of the active arcs. However, the 
Nankai Trench is much less well marked 
topographically than the typical foredeeps, 
such as the Nippon Trench. Middle West 
Japan has no intermediate earthquakes, and 
the deep earthquakes are those of the trans- 
verse belt, belonging to the Shichito arc 
structure. The volcanoes are extinct. 

Analogy between Kyushu and the North 
Island of New Zealand would logically lead 
to comparing the outer zone of West Japan 
with the major part of the South Island. The 
principal belt of shallow earthquakes in New 
Zealand, associated with conspicuous surface 
faulting, follows the coast of the North Island; 
offshore is the southern end of the Kermadec 
Trench, belonging to the intersecting arc 
structure. The earthquake belt crosses Cook 
Strait and passes down the South Island near 
its center line, with gradually decreasing 
seismicity. On the Pacific side of this line 
is the depression, generally believed to be a 
geosyncline, underlying the Canterbury Plain; 
this depression is consequently the analogue 
of the Nankai Trench. On the opposite side 
of the earthquake belt are the New Zealand 
Alps, consisting largely of Mesozoic rocks of 
geosynclinal type; they are cut marginally 
by the Alpine Fault, beyond which Paleozoic 
rocks are found underlying the Tertiary and 
Cretaceous sedimentary blanket. This western 
and north-western part of the South Island 
might thus be compared with a portion of 
the inner zone of West Japan; the Alpine 
Fault is in approximately the same situation 
as the Median Tectonic Line. However, the 
Alpine Fault, like the San Andreas Fault of 
California, is primarily a strike-slip feature; 
in this respect it contrasts sharply with the 
Median Tectonic Line, which is described as 
primarily due to thrusting. 

To complete the comparison, the wide 
eastern part of the South Island, including 
the Canterbury Plain, would correspond to an 
area south from the Nankai Trench; although 
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far off shore, this area has been suspected 
by several authors of being continental in 
character. The volcanoes of Banks Peninsula 
on the New Zealand coast are closely similar 
to typically Pacific structures like those of 
Hawaii, Tahiti, and Samoa; they probably 
have no close correlatives in the Japanese 


area. 
The two principal orogenies of New Zealand, 


the post-Hokonui (Mesozoic) and the Kaikoura 
(late Tertiary) are comparable with those of 
California and of many other regions, parti- 
cularly in that the Mesozoic orogeny has left 
many traces of thrusting and folding in the 
rocks, while the younger orogeny is chiefly 
responsible for the present block structure. 
While a similar remark applies to West 
Japan, it will not apply to East Japan, or in- 
deed to any of the more active arc structures 
of the world. 

The area of block faulting which includes 
California is larger than West Japan. Point- 
by-point comparisons between such widely 
separated areas are hazardous. California 
has no offshore trench, nor any clear evidence 
of the other arc features. However, the struc- 
tures and seismicity of the Japanese outer 
zone are perhaps comparable with those of 
the California Coast Ranges. The Paleozoic 
and intrusive rocks of the Japanese inner zone 
are not greatly different from those of the 
Owens Valley region and the Sierra Nevada. 
The extinct volcances aligned parallel to the 
Inland Sea suggest comparison with the Plei- 
stocene volcanism of the eastern Sierra Neva- 
da, Owens Valley, and the Mono Basin. 

The writer was at first inclined to cor- 
relate the Inland Sea with the alluviated geo- 
syncline of the San Joaquin Valley; but it is 
at least as reasonable to suppose that the 
Japanese correlative of the San Joaquin Valley 
was compressed out of existence by the thrust- 
ing which established the Median Tectonic 
Line. Whereas the San Joaquin depression is a 
downwarp, the Inland Sea is a faulted depres- 
sion, more strictly comparable with the gra- 
ben of Owens Valley. 

The Great Basin extends from the Sierra 
Nevada east to the Rocky Mountain structures, 
which are analogous to those of the mainland 


of Asia; but it would be rash to suppose that 
the Japan Sea covers a structure like that 
of the Great Basin. 


There is no known feature in Japan closely 
resembling the San Andreas Fault of Cali- 
fornia, or the Alpine Fault of New Zealand. 
Both are primarily strike-slip features, with 
strikes generally along the major structural 
trends, which they occasionally cross at low 
angles. Those faults in Japan where there is 
predominant strike slip often are associated 
with topography like that found along the 
San Andreas fault; but their attitude toward 
the major trends, or their shorter extent, 
differentiate them from it. 

In West Japan block faulting breaking the 
surface has chiefly been observed in the struc- 
tural belt of the Akiyoshi orogeny, along the 
Japan Sea coast. This, like similar block 
faulting in other old structures in many parts 
of the world, is most readily interpreted as 
tectonic rejuvenation. Regional stresses of 
comparatively late geologic age are being 
relieved along old lines of weakness which 
were originally developed under a different 
stress distribution. This circumstance pro- 
bably accounts in part for the tendency to 
complex surface faulting, with the develop- 
ment of multiple fault traces in a single 
event, here and elsewhere in the general region 
of Japan. Two distinct traces of nearly equal 
importance were developed in each of the 
earthquakes of 1896, 1927, and 1943 (and in 
Formosa, in 1906 and 1935). 

The eastern segment of West Japan is 
anomalous. From Ise Bay eastward, the 
Median Line curves northward in approach- 
ing the Fossa Magna; the structures of the 
Outer Zone appear drawn out of position ac- 
cordingly. This local distortion suggests that 
the crust is weak there, and the results of 
explosion seismology indicate a special struc- 
ture. The western limit of the disturbed 
segment is in the vicinity of the faulting as- 
sociated with the Mino-Owari earthquake of 
1891. Geological maps suggest the Median 
Line is offset across Ise Bay. 

The deviation of the Median Line as it ap- 
proaches and crosses the Fossa Magna is one 
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of the major problems of field geology in 
Japan. Discussion revolves around points 
similar to those raised in California with re- 
ference of the deflection and possible breaking 
up of the San Andreas Fault, descending from 
the northwest to intersect the east-west struc- 
tures of the Transverse Ranges (for details 
see the writer’s ‘‘ Elementary Seismology ”’). 
Another similar problem exists in New Zea- 
land; the Alpine Fault is deflected eastward 
on approaching Cook Strait, and diverges into 
branches which have been variously correlat- 
ed with the faults of the eastern North Island. 

The zone of deflected trends is sharply 
limited on the east by the Fossa Magna, 
beyond which we should speak of fragmenta- 
tion rather than of deflection. The western 
limit is in the narrowest part of Honshu. 
The faulting of 1891 is somewhat east of this 
limit; on the surface it is not a single frac- 
ture, but a complex of fractures with varying 
trends. In the field, the writer saw topogra- 
phic features, west of the faulting of 1891, 
which strongly suggest geologically young 
displacements along known faults. 

The presence, in the western part of the 
disturbed segment, of the epicenters of the 
transverse belt of deep-focus earthquakes, 
indicates a major transverse fracture belt, 
extending to great depth, but reaching the 
~ surface in a complicated and perhaps shifting 
series of breaks. This is essentially the in- 
terpretation of Exara, in his paper on the 
geotectonics of central Honshu. With a wealth 
of geological detail, he represents the facts 
in terms of a relative northward displacement 
of the Shichito arc structures, with fracture 
at the surface indicated by the Fossa Magna, 
dipping westward as indicated by the belt of 
deep earthquakes. 

The Fossa Magna is a structural boundary 
not obviously like any other known feature. 
There is every evidence of considerable dif- 
ference in tectonic and geophysical conditions 
on opposite sides. However, this does not 
appear to be tectonic ‘‘ contrast ’’’ in the sense 
used by Soviet geophysicists, since there is 
no clear indication of a sharp difference in 
the rate of displacements now going on. This 
agrees with the fact that the Fossa Magna 


is not a line of seismic activity. 


§5. Northeast Japan 


Northeast Japan is perhaps the most favora- 
ble region in the world for studying a Pacific- 
type active arc. It represents a different 
stage of development from that of West 
Japan. The distinguishing features of North- 
east Japan are: 

(1) Three-dimensional structure. Instead 
of block structure with vertical faults, we 
find the typical active surface of a Pacific 
arc, dipping trom the offshore troughs under 
the land, and continuing downward and north- 
westward into the belts of deep-focus earth- 
quakes. In reality this is not a thin surface, 
but an irregular thick sheet, of probably 
complicated structure. [iiss OEE active 
throughout, but (just as in block structures) 
seismicity is confined to limited regions and 
to belts traversing it. The dip is partly res- 
ponsible for the circumstance noted by Tsusot, 
that epicenters, even of the shallower earth- 
quakes, do not align so clearly here as in 
West Japan, but show areal rather than linear 
distribution. Gravity anomalies, active vol- 
canoes, and epicenters of intermediate earth- 
quakes, are aligned longitudinally as in other 
arcs. 

(1) In contrast with the relatively simple 
alignment of the volcanoes and the other 
principal arc features is the chaotic character 
of the surface geology. Studies by KopayasuHI 
and others indicate that this was once a region 
continuous with and much like West Japan, 
crossed by the Median Tectonic Line; but in 
later geological time it has been displaced 
relatively east or southeast, breaking into 
blocks which did not all move simultaneously 
nor in the same direction. Consequently, the 
line between areas formerly on the inner and 
outer sides of the Median Line may be 
drawn differently if different criteria are used, 
particularly, if the data are of different age. 

(3) Along the Japan Sea coast is a strip now 
associated with block faulting, which broke 
the surface in 1847 (Zenkoji earthquake), 1896 
(Riku-U earthquake), and probably in 1894 
(Shonai earthquake). The existence of this 
active tectonic belt was pointed out by Omori. 
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Shallow earthquakes and block faulting occur 
in the interior of some other active arcs, in 
both Pacific and Alpide belts. 

(4) A difference in mechanism of strain 
accumulation and release is indicated by 
Tsusoi’s finding that the proportion of small 
to large earthquakes is higher here than in 
West Japan. This may be interpreted as due 
to breaking into smaller units as in other 
regions where small earthquakes and earth- 
quake swarms are relatively common. 

(5) The high seismicity off the Pacific coast 
is cut by a belt of low activity, very clear in 
Tsusol’s mapping, which includes the epicenter 
of the great Sanriku earthquake of 1933. 
Probably this belt follows a resistant trans- 
verse structure, composed of relatively unfrac- 
tured rocks which yield only occasionally, 
under strain sufficient to produce a great 
earthquake and numerous aftershocks. The 
block of the Kitakami Mountains appears to 
represent this structure on land. 

(6) The Kanto district, extending west to 
the Fossa Magna, occupies the area of inter- 
section of two arcs: the principal Honshu arc 
and the Shichito arc. The former accounts 
for relatively high seismicity in the Kanto 
district. In accordance with the three-dimen- 
sional arc structure, most of the numerous 
epicenters on land correspond to foci at depths 
of 40km. or more. Hence, even when a great 
earthquake like that of 1923 occurs there, no 
evident main fracture reaches the surface. 
The many small faults formed in 1923 were 
incidental to the general surface distortion 
of the Kanto region, which was shown clearly 
by geodetic observations. This high degree 
of surface complexity recalls that of the great 
Indian earthquake of 1897, when faulting and 
complex warping of the surface took place 
over a wide area. 


§6. Hokkaido 


The intersection of two arcs in Hokkaido 
offers a promising field for geophysical in- 
vestigation of crustal structure, provided that 
the difficulties of working in that terrain can 
be successfully overcome. There are interest- 
ing analogies with the North Island of New 
Zealand, even in geographical outline. The 


central tectonic division has suggestive simi- 
larity to the Auckland Peninsula. However, 
the active volcanism of southwest Hokkaido 
does not compare well with the Pleistocene 
volcano, Mt. Egmont, of New Zealand; and 
there is in Hokkaido no clear correlative of 
the coastal mountains of the North Island, 
with their active block faulting. In brief, the 
North Island is a region of intersection of 
block and arc structures, rather than of two 
arcs; it resembles Kyushu rather than Hok- 
kaido. 


§7. Branching 


The two principal western branches of the 
circum-Pacific active belt diverge in central 
Honshu; they approach closely again between 
northern Celebes and Halmahera, then diverge 
again. West Japan may be assigned to the 
western branch; lying between the Honshu and 
Ryukyu arcs, its position is like that of the 
great region of dominant block faulting which 
includes California, between the Alaskan arc 
and the Mexican arc. The oceanic area 
between the two branches is of special in- 
terest; at the north it lies between the Shi- 
chito and Ryukyu arcs, although the deep arc 
structures of the former, extending into the 
transverse deep-focus belt, pass under it. 
This area includes the Nankai Trench; con- 
tinental structure probably extends much fur- 
ther south, although the Mohorovicié discon- 
tinuity may be found to stand near the oceanic 
level. 


§8. Comparison with the East Indies 


If Indonesia could be investigated thoroughly 
by present geophysical and seismic methods, 
the results would probably go far toward 
clarifying the nature of both block and arc 
tectonics. The results reached by past investi- 
gators there must be considered carefully for 
correlation with work in Japan and in all 
active regions. 

The belts of negative gravity anomaly, dis- 
covered by Vening Mertnesz, have commonly 
been held to indicate the location of contem- 
porary active orogeny. The islands (such as 
Nias, Timor, Ceram) included in or traversed 
by these belts were involved in intense folding- 
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during the Miocene. Presumably the zone 
weakened at that time is now yielding under 
the general regional stress; this agrees with 
the fact that the Meinesz belts include the 
epicenters of many shallow earthquakes, some 
of them large. Shallow earthquakes else- 
where in Indonesia are mostly associated with 
block tectonics. 

The principal Meinesz belt follows the 
Sunda arc, which is the least complex struc- 
ture of the region; the arc extends from 
Sumatra and Java along the chain of the 
Sunda islands and round the Banda Sea. 
Broadly speaking, it is a typical active arc, 
with all the characteristic features; but some 
of these features are better developed in one 
part of the arc, and some in another. In 
certain respects the Sunda arc compares better 
with the arc of the Ryukyu Islands than 
with that of Honshu. Thus, the islands of 
the Sunda arc form a double chain; the inner 
arc, including Sumatra, Java, and smaller 
islands, is actively volcanic, but the islands 
of the outer arc, in the Meinesz belt, are non- 
volcanic. The Ryukyu arc has a similar double 
chain, but off the Pacific coast of Honshu 
there are no small islands in the belt of 
gravity anomalies. 

The Honshu arc, in fact, is more closely 
comparable with the Java sector of the Sunda 
‘arc. Off Java there are no small islands in 
the Meinesz belt (although there is a subma- 
rine ridge); and the Java Trench, though not 
as significant as the Nippon Trench, is the 
principal foredeep feature of the Sunda arc. 
Large shallow earthquakes are not frequent 
in the Java sector; but active volcanoes, and 
epicenters of intermediate earthquakes, are 
numerous, and the corresponding deep earth- 
quakes originate north of Java. 

The relation between Sumatra and Java is 
somewhat like that between West and North- 
east Japan; however, Sumatra is more com- 
parable to an arc structure like that of the 
Ryukyu Islands. Sumatra has some definite 
characteristics of block tectonics, such as the 
high scarps along the southern coast, and the 
longitudinal internal rift, which the displace- 
ments in the earthquake of 1892 indicate is 
a right-hand strike-slip feature. (However, 


this rift is largely filled with volcanic pro- 
ducts, and resembles the rifts in the North 
Island of New Zealand and in Kyushu). 

North of the Sunda arc there is an important 
tectonic division, across a line passing through 
central Celebes. To the west and northwest 
is a large stable area, with relatively few 
earthquakes, except near the coasts of the 
Strait of Macassar, which is an internal frac- 
ture in the stable mass. All the earthquakes 
on this side of the division are shallow, ex- 
cept for deep earthquakes belonging to the 
Sunda arc. On the other side is the terminal 
section of the Sunda arc, in the Banda Sea; 
and another active arc descends from the 
northeast along the volcanic eastern peninsula 
of Celebes. Celebes, with its environment, 
presents a degree of complexity in three di- 
mensions comparable with that of Honshu and 
its environment. Interpretation in this ex- 
tremely interesting area is rendered particu- 
lary uncertain by the lack of precision att- 
ainable in the location of epicenters, as well 
as the absence of any direct geophysical data 
on the depth of the Mohorovicic discontinuity 
or on other details of crustal structure. 


$9. West Indies 


The Caribbean active arc has been much 
studied by geophysical means. Epicenter 
locations in this region have become reasona- 
bly accurate of late years; but general seis- 
micity is only moderate, so that most of the 
important large earthquakes are known only 
from non-instrumental and historical records. 
The eastern front of the arc is double, with 
an outer non-volcanic island (Barbados) and 
an inner volcanic arc. 

Negative gravity anomalies follow the outer 
arc, beyond which there is no deep trench. 
The northern limb of the arc is quite different; 
there is a very deep external trench, but there 
is no active volcanism, and the earthquakes 
are nearly all shallow. Both northern and 
southern limbs are classifiable as of block 
structure, and both are cut by important 
longitudinal strike-slip faults. 

While the eastern front of the Caribbean 
arc may be compared with the Ryukyu and 
Shichito arcs, the remainder of the structures 
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differ so much from those of Japan that 
generalizations about arc tectonics derived 
from observations in the West Indies can be 
applied in Japan only with great caution. 


§10. Himalayan arc 


Northern India is potentially one of the 
most important regions in the world for the 
study of arc tectonics. The accessible geology 
has been well investigated, and gravity obser- 
vations are plentiful. There is still much 
need for studies of crustal structure by the 
various seismic methods. 

The Himalayan arc strongly resembles 
those of the Pacific system, except for its 
generally lower activity; in this it is typical 
of the system of Alpide arcs of Europe and 
Asia. The Ganges alluviated depression is a 
foredeep; along its northern margin is a belt 
of negative gravity anomalies, in which are 
located the epicenters of large shallow earth- 
quakes. The proportion of large to small 
shocks is relatively high, which suggests com- 
parison with regions of block structure. 
Volcanism is slight, intermediate earthquakes 
associated with the arc proper are few, and 
no correlated deep shocks are known. 

The remarkable series of intermediate earth- 
quakes under the Hindu Kush are not evident- 
ly related to the Himalayan arc. At the other 
end of the arc, the great earthquakes of 
Assam and Tibet are related to complicated 
structures which are still imperfectly known. 

Nothing associated with the Himalayan arc 
suggests close correlation with the tectonics 
of Japan. 


§11. Italian arc 


Of the arcs of the Alpide system, the Italian 
arc most clearly shows relation to those of 
the Pacific. In ‘‘ Elementary Seismology ’’ it 
was pointed out that this parallelism is espe- 
cially evident when referred to double arcs, 
like the Ryukyu and Sunda arcs. Epicenters 
of shallow earthquakes follow the chiefly non- 
volcanic outer arc of the Apennines; the 
well-known active volcanoes of Italy are on 
an inner arc. Intermediate and deep earth- 
quakes are known down to depths of 300 
and even 450 kilometers. The depression of 


the Tyrrhenian Sea is correlative with such 
interior depressions as the Weber Deep in the 
Banda Sea interior to the Sunda arc. 


§12. Correlation with tectonics of the Soviet 
Union 

In section 2 reference was made to the 
classification of tectonics in terms of geosyn- 
clines and platforms. This classification has 
been applied in detail, and with excellent 
results, to the tectonic analysis of those parts 
of the Soviet Union which have been selected 
for intensive geophysical investigation in the 
past: the Crimea, the Caucasus, Turkmenia, 
and Central Asia (especially the region ex- 
tending from the Tian Shan to the Pamir). 
The Soviet Union includes typical Pacific 
active arc structures, in Kamchatka and the 
Kurile Islands. These were subjected to in- 
tensive exploration during the IGY program, 
but results are only now beginning to be 
published, so that discussion is premature; 
the opportunity for comparison with Japan 
may be awaited with great interest. 

Platform areas in the strict sense are lack- 
ing in Japan; they may be observed on the 
adjacent mainland of Asia, even though some 
deep earthquakes related to the Pacific active 
arcs originate under the platforms. The block 
and arc structures of Japan may be regarded 
as ancient platform material in different stages 
of being broken up and transformed by geo- 
synclinal processes. 

The geosynclinal zone along the south- 
western border of the USSR is on the northern 
front of the Alpide structural belt. In the 
Black Sea, under the Caspian Sea, and in 
Turkmenia, there are tectonic downwarps, 
in which the Mohorovicic discontinuity 
descends to 40-50km below the surface, and 
there are local areas from which the upper 
or ‘‘granitic’’ layer appears to be absent. 
Under the Caucasus the Mohoroviéic discon- 
tinity, generally near 40km deep, descends 
locally to a depth of 60 km. 

In most of these areas typical arc features, 
and block tectonics fracturing the surface 
along active faults, are both absent. However, 
in Turkmenia the Transcaspian depression is 
in the position of a foredeep north of the 
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Kopet Dagh arc; and the epicenters of the 
Ashkhabad earthquake series were located 
between the mountains and the foredeep as 
in other arcs. No large gravity anomalies, 
no active volcanoes, and no intermediate or 
deep earthquakes are associated with this arc. 

The area extending from the Tian Shan 
across the Garm district to the Pamir is at 
the southwestern end of the long Pamir- 
Baikal seismic belt, the trend of which di- 
verges here from that of the series of Alpide 
arcs. The southern Tian Shan, especially, is 
a former platform region which in compara- 
tively late geological time has been broken 
into blocks by the stresses associated with the 
Alpide geosynclinal movements of the Pamir. 
(Some recent evidence even suggests consi- 
derable strike-slip displacements in addition 
to vertical movements, in the Tian Shan). 
The tectonic character of the area is not 
unlike that of West Japan. Both northern 
Tian Shan ranges have ‘‘ roots’’ represented 
by depression in the Moho; under the Pamir 
both the Moho and the intermediate (‘‘Conrad’’) 
discontinuity descend. 

The Pamir region presents some typical 
arc features, although volcanism is minor. 
However, there are large negative gravity 
anomalies. Intermediate earthquakes at depths 
of 100-150km occur under the southern Tian 
Shan; if referred to the Pamir these seem 
anomalously placed. As already remarked, 
the numerous intermediate earthquakes under 
the Hindu Kush present a very special pro- 
_ blem. 

The Garm district, between the Tian Shan 
and the Pamir, is a region of relatively fre- 
quent small earthquakes at quite shallow 
depth. These are apparently connected with 
a local small-scale block structure ; volcanism 
is not involved. 

In ‘‘ Seismicity of the Earth’’ it was sug- 
gested that the entire Pamir-Baikal seismic 
belt is at present subject primarily to block 
tectonics. This corresponds with the statisti- 
cal result that the proportion of large to small 
shocks appears to be higher there than in the 
Pacific arcs. 


§13. Concluding note 


The statistical approach, noted in the pre- 
ceding paragraph, is likely in future to become 
a valuable means of distinguishing between 
areas of different tectonic type. Tsusor has 
already applied it with success to the distinc- 
tion between West Japan and Northeast Japan, 
and to further more detailed subdivision. 
Naturally the method can be applied only in 
areas of high seismicity, or where there are 
plentiful historical data extending over many 
centuries. 

The critical reader will have noted that 
there is no sharp distinction between block 
and arc regions, but rather a gradation. At 
one extreme are typical Pacific arcs in which 
the evidence of block faulting is lacking; at 
the other are regions like that of California 
in which most of the arc features are absent, 
and the block tectonics is characterised by 
large strike-slip faulting which tends to destroy 
the former alignments. Northeast Japan and 
West Japan are not quite at these two ex- 
tremes; but their opposite characteristics are 
clear enough so that useful future results 
may be expected by comparing their features 
with those of arc and block tectonics respecti- 
vely. 
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Ramsey’s hypothesis on the origin of the Earth’s core is applied to the 
model that the D-layer and the core are composed of the molecular and the 
metallic phase of MgO, respectively. The approximation of tight-binding 
method is employed to investigate the electronic band structure of MgO crystal. 
It is resulted that the pressure-induced transition of MgO to metallic phase 
occurs at a pressure of about 1.2 10° bars. 


1. Introduction 


The suggestion that the large increase of density at the boundary of the 
Earth’s core is due to a pressure-induced transition from a molecular to a 
metallic phase, was first proposed by Ramsey (1948). He abondons the iron- 
core hypothesis and assumes that the mantle and the core have the same 
chemical composition, for example, such as olivine. Olivine, however, may 
not seem to compose the D-layer, since the 6 (=K/d, where K and d denote 
incompressibility and density, respectively) of olivine is not compatible with 
that of the D-layer deduced from seismic data. 

Recently Marnunkunt (1953) and Shimazu (1959) suggested the process 
MeSiOs 2 2MgO+SiO:2 within the C-layer, which demands the D-layer com- 
posed mainly of MgO. In fact, it will be shown in the present paper that 
the elastic properties of MgO agree faily well with those of the D-layer de- 
duced from seismic data. 

Originally, Ramsey put forward his hypothesis to account for the density 
of the terrestrial planets. On his hypothesis the pressure at the boundary of 
the core should be characteristic of the chemical composition —— MgO in 
our model. Although it is important to determine whether the dimensions 
and masses of the terrestrial planets can or not show that this hypothesis is 


tenable, the possibility of the metallic transition of the constituent material at 


3 


the boundary of the Earth’s core gives a determinable evidence to his hypo- 
thesis. 

The critical pressure ai which a meiallic transition occurs depends upon 
the variation of the forbidden band (ihe energy gap between the conduction 
and the valence band) wiih pressure. i is desired thai ihe elecironic band 
structure of MgO crystal under a high pressure of about 10° bars is in- 
vestigaied. Since at pressni we can hardly expect to carry out it experi- 
mentally, although the recent development of shock-wave technique shoud be 
noted (Alder and Christian (1956)), the theoreiical investigation is the only 


way to do so. 


2. The electronic band structure of MgO crystal 


MgO crysial consisis of two kinds of ions, Mg** and 02>, which have 
closed shell structures. We suppose that the uppermost valence bands of 
MgO crystal have 6N siaies, and ihey are completely filled with 6N electrons, 
since ihe crysial is an insulator. Then we may assume that the wave func- 
tions of these bands are approximately given by LCAO (linear combination 
of atomic orbitals) Bloch functions formed from suitable atomic orbitals which 
have the characier of 2p-orbitals of O27 ion. 

The foundamental problem in the determinaiion of ihe electronic band 
structure of a crysial in the ordinary Hartree-Fock approximation is to find 
the eigenvalues of the equat'on. 

Fop(x, k) = €p(k) da (x, k) () 
(p(x, k) are the one-electron wave functions (x stands for both space and 
spin coordinates), k is a vecior in the reciprocal laiilce (body-centered cubic 
in the case of MgQ) and the index ~ denotes the Brillouin zone number in 
k-space. op is the effective one-electron Hamiltonian operator in Fock’s 
approximation. To obiain an approximaie solution of Eq. (1), we suppose 
that the ¢»(x,k) may be expressed as LCAO. For the valence band of MgO 
erysial containing N molecules, ihe atomic orbitals selecied are the 12 free 
ion spin-orbitals consisting of 6N 2p spin-orbiials on the O?2- ions and 6N 2p 
spin-orbitals on the Me** ions. These orbitals will be denoted by ¢p(x, @), 
where the vecior @ specifies ihe position of the ions, and ihe index fp goes 
over all spin-obitals associated with the ion at a. The radial paris of ithe 


orbitals used in present calculation are given by 


P, = 0.07800 7 exp(— 0.714 7) 
+8.51793 r exp(—3.412 +) + 1.6649 r Exp( ise 7) (2) 


for the 2p orbitals of O2- (Watson (195%)) and 


IE = MoO % exp( DAS i) 
+22.795 7 exp(4.808 7) +21.348 7 exp(—7.9907 7) (3) 


for the 2p orbitals of Mg** (Léwdin (1953)), respectrively. For brevity, in 
whai follows we condence our indices (~, a) and (p, 2) inio the single sym- 
bols a and k, and write 
p=eC' (4) 
where @ and g are row mairices. C’ is taken as a square (6N X6N) mairix 
so that Eq. (2) defines a set of 6N crystal spin-orbitals. The row mairix @ 
is ordered such that ihe first 3N elemenis are the spin-orbiials on ithe O?- 
ions, and ihe last 3N elemenis are ihose on ihe Mg ions. Lei ihe indices 
vu, v and 2 denoies the orbtials on the O?- ions, and 7, m and 7 denoies the 
or bitals on the Mg? ions. The indices a, 8, 7 and o will be used to denote 
any orbitals on ithe O7- and the Mg?* ions. 
The set of orbitals @ are not all orthogonal at any finite lattice paramier. 


An orithonorml set @ are introduced by Lowdin transformaiion (for example, 


Lowdin (1956), Callaway (1958)). Thus 


O=oL (5) 
where 
L=(+S)- 
fat ¥ 
7 Sap =\ba* (x) ba (x)dx— dup | (6) 
Kq. (4) is expressed by the set @ as 
gp =0C 
C+C=CCt=1. (7) 
Using Eq. (7) in Eq. (1) our fundamental approximat’on is expressed as 
CtFC=e (8) 


where € is a diagonal matrix, and the matrix F’ has element F., defined by 


Fog =Hapt 2 Cuvt Gap? Coe, (9) 
cy 


where 


Hua =\be*(x) He (x)dx, 


| 
Gap? = \Cba% (x1) * (2) Gop (ar) bs (22) | (10) 
— ba* (1) by* (X2) Gbs (X1) bp (X2) \dxidxz ) 
where H is the Hamilionian operaior for an electron moving in the field of 
the ion-cores, and G= (1/2) n—r. 
Let B be a uniiary mairix which iransform F, ihe mairix represeniation 
of all Fug, to a new mairix €’, where 


€'= BY FB, (11) 


and Bis chosen such that the submairices €/(g,m) and €’(m, #) are null 
mairices. Then a sei ¥ of orthonormal aiomic spin-orb’ials is introduced 
according io 

A= OB. (12) 
¢ is expressed in ierms of ihe set x as 

p=xV. US) 
where V is a unitary mairix, which reduces €’ to diagonal form with the 


eigenvalues €»(k) of Eq. (1) on the diagonal. Thus 


e=Vte'V (14) 
To find the eigenvalues we first consiruct the Bloch functions 
Np(r, k) =N1? 3) xp (1, # exp Zztp-k) (15) 
: in 
Then 
Pp(X, k) = IV 12>)? g(x, k) Uy (k) (16) 
where 
Ugo (k) = N"? exp(—2atk -) Va(e, &) (17) 
Using Eq. (17) in Eq, (14), we arrive the secular equation 
det Epg (k) a Ovat (k)| = 0 (18) 
for the eigenvalues €»(k). In Eq. (18) 
Sn (k) = \ my* Ce, B) Fovtta (x, k) dx, ] a 


= Eng! (pe, (6) + Dated ee, far A)exp(2nia-k), 
A 
where 2 is a vectors which maps out the O?- ion lattice from an arbitrary 


chosen O77 ion at 4. 


The e’ is obtained from Eq. (11) by treating all non-diagonal elements 
of F as first order small quantiiies and applying perturbation theory (Grimley 
(1958)). The results can be written in terms of the matrix elements defined 


in Eq. (10) and the matrix R with elements Rag defined by 
Rip= Ba Bret (20) 


For numerical calculation, it is more convenient to express Eq. (21) in terms 
of matrix elements formed with the given non-orthoganal spin-orbitals. A 


new matrix corresponding to R is defined by 


Rapg=>Ba,Byat Bag=DiLayBye 
yv 


We can separate from R a new matrix P representing the contribution from 
the O-Mg interaction, i.e. Pus=Tu.—Ruv, where T=(1+S;)-1, and S: is 
the submatrix of S consisting of all elements of the form S,,. In terms of 
these quantities we obtain 

Eun = Hyy —2G ps?* Pay = 22 Gap Pym 2G? Pane 


m m 


in SIP Pa a8; 1 Ela) 4= i Senor 1 >a Sy itive (22) 
x oN 


where Hag and Gag" stand matrix elements like Hzg and Gag?’ in Eq. (10) 
but with the set g replacing the set ®. 
We can express all non-vanishing elements €pq/ (ge, +) as fallows: 
Cla (pe, fe) = €! yy (ee, ft) =’ ee (pt, pt) = Ao 
C! ne (pe, pet A) =C! yy (pe, eta) = Ai a= (7, r, 0) 


lee (ft, et a) = Ay! A= (7, 7, 0) 23) 
Cle (ee, uta) Sean tt A= Ge 1 Q) ( : 
lve (pe, peta) =As AS (oF Q, 0) 


Here 7 is the cation-anion distance. The numerical calculation are carried out 
for three values of r=2.10, 2.00, 1.80. A. These are tabulated in Table 1. 


Table 1 The values of €'pq in three cases of y (see Eq. (23)). 


€! ng y=2.10A y=2.00A y=1.80A 
Ao 0, 205605 0.179944. 0.100703 | 
A; | 0.001162 0.009200 0.010147 | 
A,! —(), 009746 —0,013314 —0 025088 
Ay | 0.017948 0.028438 0.041487 | 
As 0.001628 0.005233 


0.006604. 


Using Eq. @3) m Eq. (18) we now have 


Exe (k) = Ao+4Ai (cos 2z7Rz cos 2z7R,+ cos 27 rks cos 2x rk) 
+4A,cos 2rrk, cos 2xrk2+2Asz cos 4rrk, (24) 
ery (k) = —4A,"'sin 2rrkz sin 2xrk 


and the other elements follow from (24) by interchanging x, y and z. The 
secular equation (18) is easily solved for several prominent directions in k- 


space, 1. e., (k, O, O), (k, k, O) and (k, k, k). The final results are showed 
in Fig. le 
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Fig. 1. Band structure of MgO in three cases of y along three directions 
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To obtain the transition point, one may proceed in the following : start 
with MgO crystal in the ionic state and determine when the energy gap disap- 
pears. The maximum and the minimum energy values of the valnece 
band are plotted against 7 in Fig. 2., which shows the broadening of the 
band with y. Since the evaluation of the conduction band of MgO crystal 
can hardly be made by the approximation of tight-binding, we assume that 
the conduction band braodens within the forbidden band as much as the valence 
band. At the ordinary lattice parameter y=2.10 A, it is found by using 
optical method that the energy gap is about 9e. v.. With the assumption 


above cited, thus it is found that the energy gap disappears at 7= 1.87 A 
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Fig. 2. The electronic band structure of MgO as a function of 
Conduction band and valence band are crossing at =1.87A. 


3. An equation of state of MgO crystal 


To obtaine the pressure of MgO crystal corresponding to the transition 
point, we must express the lattice energy of MgO crystal as a function of 7. 


According to Born and Mayer’s expression, the lattice energy & is shown by 
E,= — (amz2*e?/r) 
+ (MC 22 1472-1914 M12) (Cel? + Care*2lP?ye—@'r/?) (25) 


per molecule, where aw is Madelung’s constant (1.7476 for the NaCl type :n 
MgO), M and M’ are the numbers of unlike and like neighbours of each 
ion, respectively (6 and 12 for the NaCl type), @’ is the ratio of the distance 
between like neighbours to that between unlike neighbours (2 for the NaCl 
type), 6 is arbitrary constant chosen to have the value 10-! erg, 71 and 72 
are what shall be termed the radii of the ion 1 and 2, Ci; is afactor calcula- 
ted by Pauling and is expresed as Cij;=(1+(ei/ni)+(2;/ni)] where 2: is 
the valence of z-ion, and the numbers of valences electrons in the outer 
shell of 7-ion. 

Acoording to Huggins and Mayer’s method (1933), the unknown pata- 
meters 71, 72 and ¢ in Eq. (23) are dertermined by using some experimental 
values of MgO. ‘The final resultsis expreseed as 


E,= —76.76605/ (7/1) +1385,24705 exp (— 4.8837 (7/70) 


+30660.44466 exp (—6.9066(7/7)) (26) 


By differentiating Eq. (26) once and twice we obtaine pressure ~, and in- 


compressibility Ay as the functions of (7/7). And the relation between 


Git) | p(.l0 dynes/em) K; (€.10 dynes/cm) d (gr/cm) 
| 1,00 0.000 1.64 3 
0,99 0,052 1,81 371 
0.98 0.109 il is, 3,82 
0.97 0.174 2.19 3.94 
0.96 0,245 2.42 4,08 
0.95 0.324 2.66 4.20 
0.94 0.413 2.92 4.33 
0.93 0.511 Bogall 4.48 
0.92 0,623 Soe 4 62 
0.91 0, 741 3.85 4,78 
0.90 0.887 4 24 4 94 
0.89 1,033 4.61 5.10 
0.88 1,194. 5.08 5,28 
0.87 | 1,376 5.63 0. AT 


Table 2. Variation of density d, incompressidility K, and pressuse p & MgO 


versus change (1/ro). 
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Fig. 3. Distribution of density (lower curve) and incompressibility 
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density d and (7/ro) is obtained by using the atomic weights of Mg and O. 
These values are given in Table 2. To examine the model that the D-layer 
is composed of MgO, these calculated values and Bullen’s values (1953) are 
shown in Fig. 3, where we find that the maximum discrepancy between both 
values 6 (=Kr/d) is about 20%. We must note that the accuracy of Eq. 
(26) decreases as the pressure increases (Wada (1959)) and the calculated 
values show the isothermal and, on the other hand, Bullen’s values the 
adiabatic. At any rate, we may conclude that Eq. (26) will be enough good 
approximation to estimate the transition point, since the process in evaluating 
the electronic band structure is not so accurate that the discrepancy above 


cited playes a heavy roll. 


4. Pressure-induced transition of MgO crystal 


to metallic phase 


From sec. 2 and sec. 3 we find that the pressure corresponding to the 
transition point (7= 1.87 A) is about 1.2<10® bars, which is the pressure of 
the depth of about 2600km in the interior of the Earth, as shown in Bullen’s 
table. The result obtained, in the other word, shows that the energy gap 
between the conduction band and the valence band of MgO crystal disappears 
at r=1.87 A. In this situation the valence electrons in ordinary ionic state 
are made free, and the electronic band structure becomes metallic. We 
must note that this transition point is not such as suggested by Ramsey. In 
the former no density-jump occurs, while in the latter density increases 
discontinuously. To estimate the latter type of transition, we must obtain 
the lattice energy curve of the metallic phase. Of course we may expect that 
the minimum point of the energy curve ot ihe metallic phase could be near 
the transition point given in the present calculation. It is, however, not 
necessary that both transition occur at asame pressure. This case that the 
transition with density-jump at a high pressure than the transition without 
density-jump, is very interesting. In this case the electronic conductivity 
increases abruptly at an certain depth in the deeper part of the D-layer and 
the density increases discontinuously at the boundary of the Earth’s core. 
The suggestion that the electric conductivity will increase abruptly in the 
case of the transttion to metallic phase, is made by Mott (1956). Mott’s 


opinion is as follows : if an electron is removed from valence band to 
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conduction band, a mobile positive hole is formed. The minimum energy 
necessary to form a free electron and positive hole will decraase as the 
lattice parameter comes closer, but cannot tend to zero. since an electron and 
a positive hole attract each other with aforee of which the potential energy 
for large y is of Coluomb’s type. ‘Thus small number of free carriers of 
positive sign is impossible in the ground state, since the free earriers of 
positive and negative sign are bounded to form pairs. On the other hand, 
when there is a large number of free carriers any pair of charged particles 
may be expected to attract each other with a force derivative from a screened 
potential, which does not necessary lead to a bound state. Thus the transi- 


tion to a state showing metallic conductivity will be sharp. 


5. Conclusion 


It is found that the valence band of MgO crystal broadens rapidly as 
lattice parameter decreases, and the energy gap will disappear at 7=1.87 A 
corresponding to pressure of about 1.2 10° bars which is deduced from an 
equation of state of MgO. We may expect that the pressure-induced transtion 
of MgO crystal to metallic phase could occur at the boundary of the Earth’s 
core. A new model that the electric conductivity increases abruptly in the 
deeper part of the D-layer and the transition from molecular (or ionic) phase 
to metallic one occurs at the boundary of the Earth’s core with a large density- 


jump, is also suggested. 
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Abstract 


Variation of the ultrasonic velocities of granite in the process of defor- 
mation and fracture caused by the axial compressional stresses was observed 
under high confining pressure up to 5,000 atmospheres. At the initial stage 
of deformation, the velocities increased with the overlapped stresses, then 
reached the constant value in cases of both ihe axial direction and the 
transverse of it. In the facture range, however, the velocities of waves passing 
through in the direction transverse to the axis remarkably decreased with the 
progression of facture, while such change of velocities was scarcely observed 
in the axial direction. 

The tendency of the velocity change as siated above decayed rapidly with 
elevaiing pressure. Then, these phenomena can be explained with the assum p- 
tions of pores, closed initially and extended by pressures or stresses at the 
fracture stage, the same as the increase of strength with that of pressure and 


the extraordinary stress-strain relation reported in the previous papers. 
Introduction 


In our previous experiments” on the mechanism of deformation and frac- 
ture of igneous rocks, we researched that the rocks had the extraordinary 
values of Poisson’s ratio in the fracture range caused by the axial compres- 
sional stresses. That is, the lateral strain increases abruptly with the increase 


of the stress in this range, and the tatio of the lateral strain to the long- 
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itudinal is far beyond 0.5, suggesting that the volume of the specimen. in- 
creases with development of fracture, because the longitudinal strain varies 
almost linearly. Under elevating confining pressure these phenomena gradually 
disappear with the increase of pressure, and over two or three thousands 
atmospheres they are scarcely observed. 

On the other hand, the sirength of rocks increses rapidly with the increase of 
pressure, but the rate of the rise of the strengih has ihe same decreasing 
tendency as of ihe variation of Poisson’s ratio in the fracture range. 

These properties will be the characieristic of igneous rocks, and can be 
elucidated by a large number of pores contained in the rock specimen which 
may act like “the Griffith’s cracks’’ in case of fracturing. The volume-increase 
at fracturing can be inierpreted wiih the assumpiion ihai ihe progression of 
such pores or cracks produces ihe accumulaiion of a bulk of cavity. The 
increase of strengih wiih that of pressure can be explained by closing of ihe 
pores. 

li has been clearly explained from ihe experimenial resulis by F. Birch” 
and D. S. Hughes® thai such cavities orig*nally coniained in the rocks have 
a remarkable effect upon ihe elasiic wave velocities of them under pressure. 
The influence is stronger under a lower pressure. On the oiher hand, ii has 
been reporied*’ that the variation of elastic wave velocities of metals caused 
by the plastic deformation is within ihe limit of only one per ceni. 

Then it is expected that at the fracturing under a relatively low pressure, 
the velocities of rocks decrease appreciably, for they have also the considerable 
bulk of cavities produced in themselves. Yei, from ihe experimenis of ihe 
strain measuremenis, it is supposed that ihis decrease must be remarkable in 
the direction transverse to the applied axial stress, though the velocities of 
the axial one may scarcely change. 

F. Birch has also suggested that the factors having influence on the seism‘c 
wave velocity in the crust are pressure, temperature and the mineral composi- 
tion of it. The rise of pressure increases velocity, that of temperaiure de- 
creases it, and the increase of the composition of the basic minerals strongly 
increases the velocity in general. 

However, under the pressure wh'ch can not perfectly close the pores of 
rocks, the applied differential stresses overlapped to ihe hydrosiai‘c pressure 


will be able to affect upon the wave velocities, for they can close ihe pores 


too. 


From the above viewpoint, we studied the variation of the elastic wave 
velocities of granite in the process of deforma‘ion and fracture caused by the 
axial compressional stresses under moderaie confining pressures up to 5,000 


atmospheres. 
Experimental Methods 


Fig. 1 shows schematically the meihods of ultrasonic measurements of 
rocks under the axial compressional siress combined with high confining pres- 
sure both a) in the direction of the axis and b) transverse to the axis. The 
high pressure equipment and the full sketch of the triaxial cylinder for this 
experiment were shown in the previous paper. 


axial stress 
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P 
Fig. 1. The schematic representation of the methods for the measurements of 
the ultrasonic wave velocities under axial compressional stress combined with 


high confining pressure, a) in the axial direction, and b) in the direction 
transverse to the axis. 


Specimens used for the measurements of the axial direction shown in F ig. 
1, a) are the long columns of Kitashirakawa biotite granite 50 m/m long and 
25 m/m in diameter, and that of the transverse direction ‘n b) are the square 
pillars of the same location, 50 m/m height and 20 m/m. jn the base. 


All specimens are covered with synthetic adhesive rubber. The barium 


titanate crystals are used as the ultrasonic transducer which generates dilatational 


waves of frequency of 500 k.e.. 
Experimental Results 


In Fig. 2 and Fig. 3, the variation of the dilatational wave velocities with 
siress respectively in the axial and the transverse direction under various confin 
ing pressures is shown. In each case, the velocities show the increase with the 
increasing of the axial stress, and reach the nearly constant value. This 
will mean that the pores of the specimen have not yet been perfectly closed 
under these pressures, and thai the overlapped forces close gradually such 
residual cavities. Under higher pressure, the additional stresses scarcely affect 
on the variation of velocity, for the pores are fully closed. The ultimate 
values of velocities are nearly constant for each specimen without regard to 
the measured directions and the confining pressures. 

In the fracture stages, however, the tendency of the velocity change is 
quite different in each observed direction. The velocities of the transverse 
direction abruptly decrease with the progress of the internal failure which is 


caused by the increase of stress and with the lapse of time, though that of 


rup.| rup. 
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Fig. 2. Effect of the axial compressional stresses on the ultrasonic wave veloci- 
ties passing through the specimen in the axial direction under various confin- 


ing pressures. 
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Fig. 3. Effect of the axial compressional stresses on the ultrasonic wave veloci- 


ties passing through the specimen in the direction transverse to the axis under 


various confining pressures. 
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Fig. 4. Variation of the ultrasonic wave velocities with pressures and overlapp- 
ed differential stresses. 


the axial one are approximaiely constant up to the rupture moment. The 
change of length of specimens is negligible to this tendency, for this effect 
is in the limit of the observed error. The above phenomena will show 
that the considerable apertures are produced in the lateral direction with the 
progress of fracture. This tendency gradually decays under higher pressure, 
it is the same tendency as observed at the strain measurements. 

As expecied before, it is shown in Fig. 4 that the elastic wave velocities 
of rocks increase with the increasing of the applied axial compressional 
siresses overlapped to the confining pressure as well as with that of the 
hydrostatic pressure. The rates of the increase by the hydrostatic pressure 
and the axial stress are the same order of the magnitude, though ihe numbers 
of the observed values are insufficient to assert it. Then it may be said that 
the overlapped forces have nearly the same effect upon the wave velocity as 


hydrostatic pressure. 


Conclusion 


It will be supposed, though it may be very impudent, that the seismic 
wave passing through ihe relatively shallow region in the crust, where the 
composed rocks are highly stressed by a overlapped force and being fractured, 
decreases its velocity in the direction transverse to the applied force. In the 
deeper part of the crust, these effecis may not be observed, because the 
confining pressure is strong enough here io b2 brought no change of velocity 
by the overlapped force. 

If the above argumeni is possible on ihe earth’s crust, the variation of 
the seismic wave velocity can be observed in the progress of the shallow 
forcus earthquakes, or even of ihe somewhat deeper earthquakes, fo which the 
region exiends to the part near ihe surface, such as accompanied by the 
surface faulting. 

It will be possible to explain the velocity distribution, considering the 


stress distribution in the crust. 
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Abstract 


We investigate the variations and other properties of the micro-iremors 
at the Volcano Aso observed by the Wichert seismographs in the Volcanolo- 
gical Laboratory since 1950. In this paper, we treat the first kind and the 
second kind of the micro-tremors of the four kinds which have been clas- 
sified by K. Sassa. He studied in decade of 1930 the phenomenon thai the 
micro-tremors stop rapidly before eruptions and again increased after them. 
This time, we perceive the same phenomenon in the erupiions of Apr., 
1953, Dee., 1957, and June, 1958. Next, we can perceive that the values of 
NS-component/EW-component of amplitudes of the second kind increase 
from 0 to | at the eruptions. It is found that this disposition occured also 
in the violent eruptions of 1933. Thus, on the assumption that the crack 
parallel to the row of the present craters, which has been infered from the 
mode of vibration at the eruption-earthquake by K. Sassa, vibrates after the 
eruption, being accompanied with the magmatic reservoir, we calculate the 
azimuthal distributions of amplitudes and indicate that the transversal com- 
ponent is the same order as the longitudinal component in the direction of 


the Volcanological Laboratory. 
1. Introduction 


The jets of lava and volcanic ashes and gas, the deformation of the 
earih’s surface near the craters, the volcanic earthquake and micro-iremors, 
and the volcanic geomagnetic variations have been studied by many investi- 
gators. With respect to the Volcano Aso, many geophysical facts were re- 
vealed by K. Sassa in decade of 1930”. Particularly, he studied in detail 


the volcanic micro-tremors which appeared remarkably at the Volcano Aso, 
and pointed out ihat ihe volcanic eruption could be foretell by investigating 
the variaiions of the micro-tremors which came before the eruption. — Ii is 
already revealed ihai not only the period but also the mode of oscillation 
distinguish the first kind of micro-tremors from the second and thai from the 
detailed investigations during the comparatively silent time, the mode of 
the first kind is a Love wave type and that of the second kind is a Rayleigh 
wave iype’. This time, we will further reveal these points by investigating 
the latest volcanic actions. In the following we describe only the results of 
the analysis of recordings of ihe first kind and the second kind of micro- 
tremors by Wichert horizontal and vertical component seismographs. For, 
although there appear the micro-tremors of ihe shorter period than that of 
the above waves and they are observed by S. Yoshikawa*® near the craters, 
ihose which have the amplitudes of the same order as the first kind of 
micro-tremors at the origin damp so perfecily as to be unable io be observ- 


ed ai the Volcanological Laboratory, owing to the short period. 
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Fig. 1. Positions of Volcanological Laboratory and craters. 


2. The micro-tremors and the eruptions between Apr. 
aud June of 1953. 


The cenier of ihe eruption shified from the fourth crater to ihe first 
and the second since Sep., 1932 and afier the comparatively silent period 
between 1934~1950, we had a violent eruption of the first crater at 11732™ 
on Sep. 27th, 1953. Namely, after the small eruption on May, 1951, we 


had the silent time of two years and then the seismographs recorded ihe 
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regular waves of the first and the second kinds of the amplitudes of 0.5 
and 1.5 since Apr. 14th and on 25th they interrupted rapidly and after 
calm state of three hours, we had a violent eruption of the first crater and 
abundant quantities of lava blocks of diameter of two feet were ejected, 
some of them being throwa 600m. The amplitude of the first and the 
second kinds of micro-tremors which were recorded by Wichert seismo- 
graphs at the Volcanological Laboratory during the active period, are shown 
in Fig. 2. The arrows indicate the eruption times when the lava blocks 
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Fig. 2. Mean amplitude of micro-tremors. 


were emitted. As known from this figure, the first crater exploded also at 
Apr. 28th, 29th and May 4th. The mean amplitude of the micro-tremors 
of the first kind which was | #, and that of the second which was 1.5 yu in- 
creased 2.54 and 3~7 # after the eruption, particularly, the second kind 
remarkably increased. Such a interruption of the micro-tremor just before 
the eruption is a character of that in April, and this fact gives the power- 
full clue for the foretell. We see the same disposition also in decade of 
1930. 


3. The micro-tremors and the eruptions from 
July to Sep. of 1955. 


The micro-tremors of the first kind which have been silent since July, 
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1953 increased on July 19th and ihe second increased on 24th. The second 
which had been the regular wave gradually became the continuous and ir- 
regular wave on July 25th and the amplitude increased to 3, but, this 
time, the phenomenon that the micro-tremors interrupted just before the 
eruption was not clear and we could asertain only a irregularity in the wave 
form. Namely, we had the eruptions a few times on July 25th, particularly, 
the eruption at 1013” was considerably strong, many lava blocks being ejeci- 
ed, and then the small erupiions were repeated on July 28th and 29th. The 
amplitudes of the first and second kinds of micro-tremors in this period 
are shown in Fig. 3. 
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Fig. 3. Mean amplitude of micro-tremors. 


Also, this time, the amplitudes after the eruption were larger than that 
before it. The micro-tremors of ca. 3 continued during Sep. and then 


dampe d gradually. 


4. The micro-tremors and the euptions from Nov., 
1957 to Jan., 1958. 


After two years repose, on Nov.. 17th, 1957, the micro-tremors of the 
first kind and the second kind began to increase gradually both in frequency 
of occurence and in amplitude. The first kind interrupted on Nov. 18th 


ihe second on Nov. 28th and the eruption occured at. 19% on Dec. Ist. The 
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micro-tremors for this time .are shown in Fig. 4. 

After this eruption, while ihe amplitude of the. first kind was ca. 2 yp, 
the state of the second kind continued to rest, and the eruption occured at 
1” on Dec. 3rd, and then, the micro-iremors of the both kinds were slightly 
stronger than that before the eruption, continued in Dec. and damped gra- 


dually to April. 
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Fig. 4. Mean amplitude of micro-tremors. 


5. The micro-tremors and the eruptions during June, 1958. 


The smallness of the amplitudes of micro-tremors was a character of 
the eruption on June 24th, 1958, in comparison with the scale of the erup- 
tion, that is, the Wichert seismographs in the Volcanological Laboratory 
recorded only the micro-tremors of ca. 0.2~0.5 ~ for a few days. Howe- 
ver, the continuous micro-tremors of ca. | ~ were recorded for some days 
before the eruption by the short period seismographs of very high magnifi- 
cations in the Hondo’s Observatory near the crater and interrupted two 


hours before the eruption and increased after it. 
6. The characters of the variations of micro-tremors. 


From ihe variaiion of the amplitudes of micro-iremors of the both 


kinds which occured between 1953~1958, we can summarize the following. 


li 


1) The micro-iremors of ihe both kinds occur several ien days before 
ihe erupiion, ihe amplitudes larger than 1 y ai ihe Volcanological Labora- 
tory coniinue and ihe micro-iremors which reach a maximum in amplitude 
and frequency of occurence ofien siop rapidly several hours or several days 
before ihe erupiion. However, this disposition was not evideni only in ihe 
eruption on July, 1957. 

2) Afier the micro-iremors of ihe first kind of the amplitude of ca. 
1 w coniinue io be recorded wiihin ihe limiis of several days or several tien 
days, ihe siaie of them often reiurn io rest without ihe erupiion. For ex- 
ample, ihe end of June, 1954, the beginning of Jan., 1955, about ihe mid- 
dle of Jan., 1957, and eic.. These are distinguished from ihe dispostiion 
of the micro-iremors before the erupiion, by being accompanied wiih ihe 
second kind, and by the gradual siop of the first kind before ii. However, 


to foretell the erupiion, there are many poinis io siudy yet. 


7. The variations of the mode of vibration of 


micro-tremor of the second kind 


Irrespective of ihe aciivity of volcano, there occur the volcanic microire- 
mors which have been called ihe second kind by K. Sassa, and ihe period of 
which is 3.5~8.0sec.. The mode of oscillation was invesiigaied by means of 
ithe simulianeous observations with ihe Galiizin type seismographs, which were 


set two poinis on Oci., 1932, and which were sei six poinis on Aug., 1933. 
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Fig. 5. Direction ‘of horizontal displacements of micro-tremors of 


the second kind (after Sassa). 
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These resulis are shown in 
Fig. 5. 

Namely, ihe second kind 
seems io be a kind of longi- 
tudinal wave from the faci thai 
the direction of the horizonial 
componeni is that of the pro- 
pagaiion. The veriical com- 
pononi is the same order as 
ihe horizonial component or 
less than thai, and then they 
may be a kind of Rayleigh 
iype, but noi che ordinary one. 
Thai is, the horizonial longi- 
iudinal vibrations viriually co- 
incide in phase with the verii- 
cal, which are not in accord 
with the theory of ordinary 
These resulis 


were obiained during the silent 


Rayleigh waves. 


period when a series of erup- 
tions of 1932~1933 drew io a 
close, and the amplitude was 
small, e.i. 1 and at ihe Volca- 
nological Laboratory only ihe 


EW-componeni was recorded. 


And also in ihe case of ihe eruption of Oci., 1932, they were about ihe 
same. ‘These waves which occured iowards ihe close of volcanic aciive period 
seem io be creaied from ihe magmaiic reservoir of the considerable dimen- 
sions which is under the crater. We investigate ihe variations at erupiion 
of ithe mode of vibration of the second kind of ithe period of ca. 3~4 
sec., because ihe micro-iremors of ihe period larger than 5 sec. were 
not recorded. That is, the daily mean values of NS-component/EW-com- 
ponent of ihe second kind in Apr., 1953, July, 1955 and Dec., 1957 are 
read. The results are as follows. 


The values increase from 0.2 io 0.8~0.9 ai ihe volcanic eruptions. 


May i? , 


sec. 
Fic. 7. Portions of records obtained from the Wicherts seismographs 
fags 


on Apr., 15, 1953 and May 17, 1953 at Vol. Lab. 
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Namely, ihe NS-componeni 
which do noi occur yet before 
ihe erupiion occurs -afier ii. 
Also the mean values for the 
erupiions of 1933 are read, as 


shown in Fig. 6. As known in 


this figure lv as evident that the 15 20 25 30 5 10.15 20-25 30 5°10 15 20 
: ~ Jan. 
same phenomena iook place also mds 


in them. ‘This faci suggesis that 
a wave generates from the vibra- 
tional origin other than ihe mag- 
maiic reservoir, being accom- 
panied with the vibration of ii. 
Nexi, NS-componeni/V-compo- 


10° 15920-25250) “Sead tS 


neni observed ai ihe Volcanologi- Sep. 


cal Laboratory are ploited in Fig. 
8. We obiain the resuli corres- 
ponding to ihe case of NS-com- 
poneni/EW-componeni. Thai is, 
while the veriical componeni of 


the same order as the EW-com- 


ponent both before and afier ihe 
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Previously, K. Sassa siudied the mode of 
vibration near ihe crater of the small erup- 


/ Hondo tive earthquake on July, 1933. The primary 


shocks are shown in Fig. 9. 
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Guake of the, Pixst.Ceaten explained by the vibration of crack parallel io 


in Aug. 1933 (after Sassa). the row of the present aciive craters. 


Wa 


After ihe considerable erupiion, ihe ground will loosen near ihe reser- 
voir and then ihe above veriical crack will. become io be able io move per- 
pendicular io ihe crack plane, being accompanied with ihe vibraiion of 
reservoir. Then ihe remarkable variaiions afier the erupiion of the mode 
of vibraiion of ihe second kind seem io be due io ihe waves generaied 
from ihe above movemenis. In order that we examine if there occur ihe 
iransversal component of ihe same order as thai of the longitudinal com- 
ponent in the direciion of ihe Voleanological Laboratory, we calculate ihe 
vibrational modes of ihe waves round ihe craier generaied from ihe move- 
menis of crack of ca 1km., which coniinues io the reservoir and is parallel 


to the row of craters. 
8. The distributions of amplitude. 


There lie ihe above crack of finite dimension, a ceriain depih under 
ihe surface, however, ithe complete ireaimeni of such a problem would 
be so difficult to solve that we neglect the effect of the free surface 
and assume for simplicity that the lengih of ihe veriical direciion is 
infmiie, that is, the problem is reduced to iwo dimensional one. As we 
calculate the horizonial propagaiion of wave in the above problem, it seem 
io be valid. Using ihe ellipiic coordinate, the wave equations are wriiien 


as follows 


074 At+2 CA Od) © 
oO C2 roakiae See maa at On? a | 


Oa fe (7 
? “Ot ~ ¢2 (cosh? € —cos? 7) \ OF + aya 


x=c cosh €-cos 7 


y=C sinh ¢-sin 7 

A= hi% a + a iF 
“itl oe (Ge) tae eo} 

] 


——- = c* (cosh? § —cos? 7), 
hi? 


no) oa | 


where 
“u,v =componenis of displacement refered io curvilinear coordinaies, 


o=densiiy of isotropic solid, 


12 


4.u=Lames elastic constants. 
According to K. Sezawa,”” ihe solutions of these equations are express- 


ed as follows 


4 =3} BuHn(é, )Gn(1, em 
20 = 31 Can, 2Ga(n, ae™, | 
n=0 


where 


SOO ee OP get EE ge 
1-551 ieee an 
Gn(7,@), Gn(7, 9") are Maihiues funciions and An(€,q), Hn(&, 4’) are ihe 


soluiions of ihe following equations 


CHE? + (h®c? cosh? € — 1") Hn(, 7) =0 
d@?Hn(E, 9") 


an + Cec? cosh? € —n*) Hal, a!) =0. 


. ta 
u1, Vi the componenis of displacemenis of ihe dilaional wave and 22, 
v2 those of ihe disioriional wave are given by 


° IyBn On 


se J (€, 9) s 
Wi=— 2 pa OE Gra, De 


n=l h? /) 
a he aC é 2 aS Se Ree Sa (4) 
= 1 Vn = ! ‘Gn 3 $0 
Us = pay 2 Ane qd ) an etre 
os y Gn OHn G74 ener q' ett, | 
n=1 Rk? 0& / 


In ihis case, as the gas rich magma seem io press only veriically ihe 


surface of crack, ithe boundary conditions are given by 


normal __ ofa Oe ee ee O(1/hy) 
normal [Ad-+ 2p + Corte) fio}) = Se o 


iangeniial a ue (Ui+v2) i Oh, (a+ Uo) 


| 
siress 0& On leg =0: 
When € approaches 0, asympvo.ic relations are known such as 


Hn(E, 1 OED, _ ny Ieee 


' = = — 
Hen E89 sy wea 


ais 


Taking ihe physical conditions inio consideraiion, Mathieu funciion for 
the dilatational component is even with respect io 7, while thai for the 
disioriional component is odd. Then, ihe boundary conditions (5) are 


reduced inio ihe following expressions 


S= [3 BnGa(n, Dt = : 3 anBuGn(7, QD+> BnCy PCat We 2 2 
| 
OGn | 
ergy Ba ge +B BaCaGuCn, 2} | 
= ( _ Go UBS OGo(%, i — An 'BnOGa(”, a a On 'Bn si sin 27 
vel \ sin? 7 On 2m isin? uy On 2 “sin 7 ~ Gah, D 
as Br! Cr | Galt, 74. SD I (Cer 'Cn sin 27 0Gr(%, q') 1T'Cn O7Gnr(%, q ny | 
ncisin? 7% vt sint 7 On sin? 7 an? j 
en art (6) 
where 
et =26V Hct Cae 
an = 2 h2ce =i: E.——— hec2 T h2c2 ? 


awe Se eas poe Le 
hc? ee kc 


an 


In order io solve ithe problem approximaielly, we use Maihieu func- 
tions of Oth io 6:h orders, and negleci the higher terms in g. The Maihieu 


funciions are given by 


c@é)=1+4+ 44 cos 27, cCé,=cos 4+ 4 cos 37, 


Cé2 =cos 27+ 4( 2 cos 47 — 2), Cé3=COS Beal cos It cos 57), 
C€4s=COS 4n+a(——— cos ant = cos 61), 


: : F 2 : 
Sé,=sin 7+ 4 sin 37, Sé2=sin 27+-2- qsin 4%, 


Sé3=sin 3n+a(—sin q+ ; sin 52), Sé4=sin 4n+4( 4 sin 47+ z sin 67). 


Using these funciions, ihe normal siress and ihe iangeniial siress ai ihe 


crack plane are oxpressed by 


n=1 


Sot >) Sn cos 27 | 
| 
| 


= 1[ SuBy+ SieBot SB SiuCst SsCt tae 
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where 


+ (Sa Bo Sze Bat San Bi Sra Seca ). : 


+ (S31Bo + S32Bo+ S33Bit S3iCot+ S3sCa+ -+°**° 


te (Sir Bo + SieBot Sis Bit SysCot+ SusCyteree 


oe) 


-) 


n= 


Sn’ sin 2H 
1 


= 1{ (Su!Bo+ Sie'Bot Sie! Bet Su'Cat Sis!Cat ica 


+ (Se1'Bo+ See’Bet+ S23’Bit+ Ses’Cot+ So5'Cy+ 


Si =14 8¢ao 


Siz= — 24+ 2ae+—$- Pier ae: ar 


3 


Sis= 0 Gas+4as— 327 se ar 
Sux ~ 26, — q'B. 


Sis = -BA.-< q'B 


Sor =" — (Koi —4Aqgay 


Sea = ~ ast -5 gast47+"2 qr 


: 176 
Ses = — aut 16r+ 46 qas+ ve qr 


So4=0 
S25 =0 
S31 = 49g 
8 3 
Ss=1+—3 ga, — ar 


Pa 


Sa oo “ q+4a4+ “ qa, — 167 —->= qy 


S34= —8q'B> 
Sue= -128,—% a’, 
Sai =0 


pala | 
sin? 7 
eos 27 | 
Yeos 47 | 
ete (8) 
cos” 
sin? 7 | 
iene sin 27 | 
192 
5 
SrA IND 0, (9) 
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Sa= qd 
Sig = 1+ Gas — = qr 
24 =O 
Sis= —89'B4 
Sit! =2a0'+8¢ao' 
Six!=2as'— 3 gas! 
Sis!=2a,'— 3 qa! 
Sis! = —47'— 16 g’y’ 

3 

aA 


Sip = =a ae ax! 


Now we consider the case that uniform pressure changes periodically 
with time, and the iangeniial siress is zero. Moreover, we may pui for 


symmeiry, 


B,=B;=C,=C;=0. 
Then, we obiain ihe following simultanious equations 


Si Bo+ Si2Bo+ SisBi+ SisC2+ SisCs = So/t 
S21 Bo + S22B2+ S23Bi+ S2iCo+ S2sCs = 0 
S31 Bo + S32Bo+ S23Bit S31C2t+SasCi =0 Ie ee 10) 
Su Bot SaBrt SisBitSuGe+ SisCi =0 | 

Si1'Bot Si2'Bot Sis'Bit Sis'C2+ Sis'Cs =0. 


We obiain By, Bz, Bi, C2, Ci from these equaiions. 
When ¢ is large, ihe asympioiic expressions for Hae, @), Hn, a") 
are obiained in ihe following form. 
—ihe sinh & | 
hE, )- 
o—ike sinh § | 


Hn : Ps os 
Med) Theanh © 


} 
nok. ., t@ . = 
Thus, the radial componeni of ihe dilajtonal wave w and ihe iangeniial 


component of ihe disioriional wave Vz ai a distant are given on inseriton 


of (11) inio (4) by 
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co 18% 0H, (&, q) t 
= <2 ie /ecchte — coetey NE ee ae | 

=-» segs 8 geos(pt—hR)} 

n=0 

=-35 2 ear — — Gs. (a, q' ext 

n=1 C08 cos? 

dj sprapiin\ onsen q')sin(pt — RR) + dnsen(7, g')cos (pt — RR)}, 
2 eee (12) 
where 


Bn=an+ibn, Cn=Cnt+tdn 
R=ihe disiance from ihe origin. 
Now we put 
T=4sec. 
26/v.= 132, 
where 2c=lengih of crack. 
On resulis of calculaiiors cf dn, bn, Cn, as we obiain 
ay) = 0.02346, do = —0.01584, a= —0.001562, 
by = — 0.09339, ba=0.08791, b, =0.0008799, | Ait aaa) 
dz =0.1214, d,=0.001355, 
C2= —0.1645,  c,=—0.001651. | 


[uJ 


s The azimuihal disiribu- 
S : tion of #4, Ve in this case 
is shown in Fig. 10. 

We ihoughi previously 
that the was riched magma 
in the reservoir which was 
infered io exist under the 


craiers frequenily exploded 


in the aciive period of the 
Voleano Aso and being 


accompanied wiih ii, there 


I ; . : 
: occured the micro-tremors 
of the second kind, and 
al Seen fg ee Gila ee : 
Fig. 10. Azimuthal distributions of displace- then indicaied that ihe di- 


ments, #, and v2, at a distant point. mensions of ihe reservoir 
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corresponding io ihe period of ihe second kind is ihe order of 1 km, on the 
assumpiion ihat ihe form of ihe reservoir is a sphere®”. As a resuli of 
the above investigation, we find thai a iransversal wave propagates simul- 
taneously with a longitudinal wave during ihe shori period afier ihe consi- 
derable erupiion. The above calculations indicate that when there exisi the 
crack of the order of 1km which conneci with ihe reservoir and after ihe 
erupiton it operaies as ihe origin of waves with reservoir, ihe transversal 
wave of the same order as the longiiudinal wave propagaies in ihe direciion 
of the Volcanological Laboratory. Therefore, it is infered from ihe ap- 
pearance of ihe iransversal component only during ihe shori period afier 
the eruption thai there occur ihe rapid changes in ihe ground around ithe 
reservoir and in the distribution of ithe gas riched magma and thai they 
restore io ihe state before the erupiion wiih ihe decrease of ihe volcanic 
aciivity. As we did noi observe at many poinis for ihe long period, we can- 
not bring ithe unambiguous conclusions, however, the above mechanism of 
the wave generaiion near ihe craier can explains ihe observed resulis. 

The writer wishes to express his heariy thanks io Prof. K. Sassa for 


his insiruciions. 
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The accumulation of data of the physics of 
semi-conductors and dielectrics created new 
possibilities of studying the thermal regime 
and the thermal history of the earth. For 
instance the role of the radiative thermal con- 
ductivity was investigated (CLARK 1956, Pres- 
TON 1956) which in addition to the lattice ther- 
mal conductivity must lead to an intensive 
heat transfer in the deep interior of the earth. 

The investigation of the pattern of the lattice 
part of the thermal conductivity led to a con- 
clusion about the presence of the resultant 
thermal conductivity minimum at 50 to 
100km depths of the mantle. This is of a 
decisive importance for the study of the earth’s 
thermal regime. The low thermal conduc- 
tivity of the earth’s upper layers hinders an 
intensive loss of heat and provides for the ac- 
cumulation of heat and for heating of the 
deep interior of the earth but not its cooling. 
Due to the minimum of the thermal conduc- 
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tivity the temperature gradient has an in- 
creased value which is sufficient to explain the 
observed lowering of the velocities of seismic 
waves at these depths (Lusimova, H. A., 1958, 
1959). 

The new data about the properties of ionic 
semi-conductors and non-metallic materials in 
the range of high temperatures indicate the 
existence of one more process of energy 
transfer by the excited states of atoms or ex- 
citons. The possibility of the existence of 
exciton states in a given material is established 
by the absorption spectra of the material. 

Compounds MgO and Mg:,SiO, are the main 
components of the typical earth rocks and 
meteorites. Kincery (1954) measured the ther- 
mal conductivity of these materials at high 
temperatures. Fig. la shows that at 7T> 
1000°C the thermal conductivity passes through 
its minimum and begins to grow with the in- 
crease of J. The deviation from the lattice 


ant’ il ; 
thermal conductivity law = or increases 


7 T/A 
with the increase of temperature (Fig. 1b). 
Let us compare these data with the data con- 
cerning the absorption spectra of these materials. 

The absorption spectrum of olivine was 
studied by CuarK (1957), special attention be- 
ing paid to the existence of the transparency 
interval providing for the radiative energy 
transfer. The presence of several absorption 
peaks in the transparency interval remained 
unexplained although a hypothesis about their 
possible exciton origin was assumed by CLARK 
and Jamieson and Lawson (1958). 

Reming and Henstey (1958) found strong 
and weak absorption peaks in the spectrum 
MgO at T=—170°, +25° and +150° and as- 
cribed them solely to excition state. Exciton 
peaks are found in the spectra of a number 
of other substances, such as KI, KBr, InSb, 
Cu,0, Pblz, TeBr, Ge and others (Haken, 1959). 
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In atomic semi-conductors these lines can be 
observed only at low temperatures when the 
background absorption is weakened, whereas 
in ionic crystals the exciton absorption is so 
strong that it can be observed even at high 
temperatures (MarTIENSSEN, 1959). CLARK’S 
observations are made at room temperatures. 
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Fig. 1b 


The distinguishing of exciton lines is facili- 
tated by the fact that a definite type of exci- 
tons—non-localized excitons— gives lines narrow 
in their spectra with the frequencies satisfy- 
ing a hydrogen-like law: 
Wes 


Yn = Vo— 


(1) 


We found that several narrow lines in the 
absorption spectrum of olivine can be identified 
with the exciton state of a definite quantum 
number 7=1, 2,3 assuming that v..=29000cm-!, 
Rex=25800 cm-!. These are lines with fre- 
quencies 3200, 22600, 25200, 29100cm-!. The 
condensation of lines occurs in the direction 
of the ultra-violet absorption, the interface 
frequency of the series coincided with the edge 
frequency of this band which is equal to 
29100 cm=}, 


The value of the constant R-x=25800 yields 
a correct value for the given mass of exciton 
Rezno'm 

R 
where R=1,097-10° cm-! is the RYDBERG’s con- 
stant in wave numbers, m=9.11-10-** g is the 
electron mass, and m is the refractive index. 
The diameters of excitons can be calculated: 


M.*= ~2m 


2 


h 

du= 27 roa 

For instance for m=2 and 5, we obtain d:= 
6A and d;=150A. 

These numbers do not contradict with the 
diameter values established for other semi- 
conductors. The considerations mentioned 
above permit us to believe it plausible that in 
olivine exciton states exist with the energy of 
excitation: 

Fi=0-4e.V., Ex=2-8e.V., E3=3-12e.V. 

There is a difficulty as Dr. Ciark said to 
me. It is difficult to fix exactly the interface 
frequency. As usual the data are scanty. 
One must make new experiments in this direc- 
tion. 

In the light of the above stated the energy 
E=1.19 has no relation to excitons as was 
erroneously assumed by JAmiEson and Lawson 
(1958). It is associated with a wide absorp- 
tion band at frequency »=10,000cm-! which 
is due to the presence of Fe ions. All the 
excitons should make a contribution to the 
thermal conductivity, therefore 

Avz=A\(Ex) + A2(E2)+43(Es) ++ +++ 
where 


nh? a 


ool E(B) a] 

ICT ONS FOL : 

1=1)/T being the length of the free path of 
excitons. In normal conditions practically 
only the first number 4; with the lowest ex- 
citation energy plays the main role. With 
the increase of temperature exciton peaks 
shift in the direction of longer waves and 
lower energies (MARTIENSSEN, 1959). There- 
fore at high temperatures the contribution 
from the subsequent terms must increase and 
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the role of the exciton thermal conductivity 
must increase. The rate of the peaks’ move- 
ments with the increase of JT in the direction 
of greater energies in MgO is =, 10-6. Vo 
(Rmitinc and Henstey, 1958). 

The contrary effect on the parameters of 
the radiative thermal conductivity 4,.= 
16om0?T*/3€ is caused by temperature. The 
absorption coefficient € of the radiative energy 
depends on T and » in a complicated way. 
At different frequencies different absorption 
mechanisms exist but they all lead to an in- 
crease of absorption with the increase of tem- 
perature and pressure and the dying away of 
the radiative heat transfer at great depth. 

In the transparency interval of olivine (0.3~ 
11)“ the absorption is due to free carriers, 
excitons and split of the energy levels of the 
unfilled electron shells of iron Fet++. The 
formula for € can be written as 


exp(——4-) 
KT 
ar ee ae 


yp? 
V> maz 


ee pt erat ee 
where & is the absorption due to the scatter- 
ing processes. Fy, is the energy of the activa- 
tion of the carriers, »: is the frequency of the 
i-th exciton absorption peak, is the 
frequency of the maximum absorption by ion 
Fet++, dy is the width of the absorption band 
due to Fe**. 

Such complicated dependence of € on vy and 
T is rather difficult to be used. For the sake 
-of simplicity the RossLanp averaging formula 


can be used, then the average value & of the 
spectrum must be calculated for the given 
temperature T and the pressure P from the 
formula 


6=60-+61 


SMSET De) {-a 
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Vmax 
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é i Art Jo Eu (1—e-“)? 
where u=hv/KT, h is the PLancxk’s constant. 
Formula (2) for € includes the value Fy the 
energy of the activation of the current car- 
riers. The energy of the activation of the 
ionic electrical conductivity predominant in the 
upper layers of the earth is studied in a 
number of papers (see Review of Prof. Tozer, 


1959). The dependence of the activation ener- 
gy of the intrinsic electron conductivity pre- 
dominant in the lower layers of the mantle on 
T and P can be written as 


OR, OEy 
aP ap (3) 


where a@ is the coefficient of the thermal ex- 
pansion; K is the incompressibility, 2% is the 
width of the energy gap at normal pressure 
and temperature. The second term allows for 
the effect of the pressure on the width of the 
zone. According to the data of Baucuan, A.S., 
OE, 


and Drickanegr, H.G., (1959), aP = —2.5. 
e. V. 


UG Apnetea Bar The third term allows for the 
effect of temperature on F, due to the thermal 
expansion of the crystal and the fourth term 
allows for the effect of the interaction of elec- 
trons with the thermal vibrations of atoms 
(phonons). Irrespective of the crystal struc- 
ture the interaction of electrons with phonons 
always leads to contraction of the width of 
the energy gap and to an increase of the 
electron conductivity. The value of the coef- 
ficient 8 is B=£o/0?/3(u)? (uw is the sound velo- 
city). The results of the calculations of the 
valus EH, (P,T) and the distribution of the 
intrinsic electron conductivity o=0oo exp (—E,/ 
KT) are given in Table I, our former results 
about the distribution of temperature in the 
earth having been used. 


E,=E%+P 


—akT 


sap L, 


Table I 

- Co 
km Ge ae r Eg Theoretical 

variations _ ee. Sak — 
750 Da lOm 2920 | 2.6 2.10-2 
1500 ALO" 3960 | 2.0 Zen 
2000 UNO! ALGO) |) 125 AKO" 
2900 ae 4300 | 0.6 DAI 


We can see a Satisfactory coincidence of the 
calculated values of o for the lower part of 
the mantle with the data of the value oa ac- 
cording to McDonatp (1957), obtained by 
means of geomagnetic variations. This indi- 
cates the plausibility of our initial assump- 
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tion about the predominant role of the intrinsic 
electron conductivity in the lower layers of the 
mantle and confirms also that our former tem- 
perature calculations are in agreement with 
geomagnetic data. 

Applying the above mentioned considerations 
we calculated the absorption coefficient of 
olivine at T=1500°K and obtained that = 
500cm-!. The curve of the thermal conduc- 
tivity of olivine against temperature can be 
best explained if it is assumed that besides 
lattice (phonon) thermal conductivity the radia- 
tive transfer and exciton thermal conductivity 
are present, the absorption coefficient varying 


from €=20cm~-! at room temperatures to é= 
500 cm-! at T=1500°K and the length of the 
free path of excitons being /=10-°cm at T= 
300°K as it is the case for a number of other 
semi-conductors (DEVYATKOVA 1957). 

The theory of the lattice (phonon) part of 
the thermal conductivity of dielectrics in the 
range of high temperatures was dealt with in 
a number of papers. Most of them are based 
on the allowance for the anharmonicity of the 
thermal vibrations of the crystal lattice. The 
ultimate formulas differ somewhat from one 
another but it can be shown (Lawson 1957) 
that they are equivalent and are reduced to 
the unique expression: 

Apn=Ao/ 37° T X77 01/2 
where X is the compressibility, 7 is the 
GRUNEISEN coefficient. Special place is oc- 
cupied by the formula of PomerANcuuxK (1943) 
who besides the anharmonicity took into ac- 
count the dispersion of the thermal vibrations 
of the lattice. It cannot be reduced to ex- 
pression (4) and is written as follows: 


= 


uey'/2 

“refi (5) 
Since the temperature of the matter of the 
earth’s mantle is very high, close to the tem- 
perature of melting Tm, we can use the rela- 
tions existing at such temperatures, such as 
LINDEMANN relation. Inserting the compres- 
sibility value from the Linpemann relation 
instead of X into (4) we obtain the following 


relation between 2ph and Tm (Keyes, R. W., 
US I518))- 


A ph= Ay o/s 


AphT = BoTm3!?207/2A/® (6) 
Bo=R?/?/372E2No'/? 

is the average atomic weight, € is the ampli- 
tude of the thermal vibrations at melting point, 
No is the Avocapro number. The value Bo 
can be considered constant to a rough ap- 
proximation. The validity of this relation is 
tested by Keyes for more than thirty dielec- 
tric crystals. A similar transformation can 
be made with the Pomerancuuk formula and 
we can obtain: 


APhT §/4=Bip?Tm?/ A? . (7) 


Inserting the values typical for the mantle of 
the earth into formula (6) we shall find that 
the values of Tm in the mantle will be less 
than those estimated by R.I. Urren (1952) and 
Juarxov (1959) and given by McDonatp (1959). 
A better approximation is obtained if relation 
(7) is used which is derived from the Poms- 
RANCHUK formula (5), and increase the average 
atomic weight of the material in the lower 
part of the mantle believing it to consist of 
the orthosilicates of iron. 

We must not overestimate the accuracy of 
these semi-empirical relations. But since our 
knowledge about the behaviour of the material 
in the deep interior of the earth is very scanty 
it is sometimes very important to establish 
that various values do not contradict with one 
another. An example of it is the above con- 
sideration as well as the comparison of the 
electrical-conductivity pattern with the distri- 
bution of temperature. 

A comparative distribution of the three 
components of the thermal conductivity in the 
earth’s mantle is represented in Fig. 2. The 
temperature distribution for the present mo- 
ment was taken from our paper (1958) for the 
calculation of the curves. We see that the 
minimum of the thermal conductivity in the 
upper layers of the mantle becomes sharper 
than it was earlier without the allowance for 
the exciton component. The thermal conduc- 
tivity at the base of the mantle exceeds its 
value at the suface by almost two orders. How- 
ever this does not practically change the in- 
ferences about the process of the thermal his- 
tory of the earth made earlier without taking 


where 
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into account the exciton component on the 
assumption that the absorption coefficient of 
the radiative energy €=10cm—!. 

Let us consider now the case of the infinitely 
great thermal conductivity in the layer of the 
mantle lying below the level of the minimum 
thermal conductivity. Our former calculations 
as well as the calculations made by McDonaLp 
(1959) showed that at the depth of 100 to 700km 
the temperature was close to Tm and in this 
region conditions seem to have been favour- 
able for convection. The moment of the ap- 
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pearance of such a layer refers to the time of 
2-10° years ago. Together with Van Zin Di 
we calculated the temperature distribution and 

the heat flow which could have existed in such 
an extreme case where the thermal conduc- 
tivity is infinitely great in the indicated layer. 
The calculations were made by Lyxyanov’s 
hydrointegrator by the method described in 
our paper (1958). The infinitely great thermal 
conductivity was reached by a complete re- 
moving of the corresponding hydraulic resist- 
ances. Figs. 3 and 4 respresent the results of 
the calculations of the heat flow and temper- 
ature distributions for different moments of 
time. We can see that the earth began cool- 
ing 1.5x10° years ago. However the heat flow 
at this period exceeds its modern observed 
value by 5 to 8 times. This indicates the 


unreality of the initial assumption about the 
existence of an intensive transfer of heat by 
convection in the mantle up till now and shows 
that such convection embracing a continuous 
belt in the earth must have stopped by the 
present time. With the end of such convec- 
tion the earth’s interior must have begun 
heating anew. This calculation can be believed 
as an ultimate proof of the absence of con- 
vection with the global scale in the earth’s 
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interior. The convection could have occurred 
slower. It is essential that in this case ac- 
cording to our calculations a less intensive 
convection cannot have led the globe to the 
state of great cooling. 

We cannot deny a convective transfer of 
heat in the earth’s mantle. It was rather of 
a local character in the form of carrying out 
of heat when the sialic material was trans- 
ferred from the mantle to the surface in the 
process of the crust formation. In this case 
the heat flow could have passed through its 
maximum several times in local areas of the 
earth (Lusimova, 1960). 
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On Conditions of Magmatism Origin and Role of Volcanic Activity 
in the Thermal Regime of Earth’s Crust. 


By 
H. A. Lupimova 
Institute of Earth’s Physics, Academy of Science, USSR 


The origin of magmatism is closely con- 

nected with the general development of the 
earth, particularly with its thermal history. 
At present a number of papers (Lusrmova 
H. A., 1952-1958, Jacogs and ALLAN 1954-1956, 
McDonatp 1959) dealing with several variants 
of the earth’s thermal history in detail are 
available. Despite some differences as to the 
initial assumptions these papers lead to a 
general conclusion about the secular accumula- 
tion of the radiogenic heat inside the earth 
and the slow heating of its interior. The 
formation of the earth by means of the ac- 
cumulation of interplanetary particles is now 
considered as the most probable initial state 
of the earth. In this case as it was first shown 
in papers of O. Y. Scumipt and his followers 
(1955-1958) the heating of the planet in the 
process of its formation could not cause its 
entire melting. 
_ The use of data of physics of solids casts 
new light on the thermal state of the earth 
and permits in particular to study the thermal 
conductivity distribution in the earth’s mantle. 
It was established that in the upper mantle 
__ (50 to 100km) the thermal conductivity decreases 
to its minimum and at great depths it rapidly 
increases with the increase of temperature 
(Lusimova 1958). This minimum of the ther- 
mal conductivity hindered the intensive loss 
of heat through the surface. Our calculations 
showed that the conditions for the cooling of 
the whole earth and for the contraction of the 
earth did not exist. 

The greatest temperature gradient was in 
the upper mantle approximately up to 100km 
deep. This caused the temperature at the 
depth from 50-200 to 500-700 km to amount 
to the temperature of melting at 2 to 3 billion 
years after the formation of the earth (Fig. 


1). In the upper part of the earth’s mantle a belt 
was formed within which the conditions were 
favourable for magmatism to originate and 
the matter started differentiation (Luprmova 
1956, 1958). The age of magmatism should 
be counted off since that moment. Before 
3-10° years ago magmatism did not exist. 

Approximately at the same time inner parts 
of the earth’s core at the boundary of which 
the transformations into the metallic state 
takes place must have melted. The recent 
computations of McDonatp (1959) made by an 
electronic computer confirm in the main these 
conclusions. 

The melting belt was not continuous. With- 
in the accuracy of the available data we can 
speak about the approximation of the tem- 
perature curve to the melting curve. The 
excess over the melting temperature is small. 
Moreover taking the heat losses due to the 
latent heat into consideration we can speak 
only about the local melting pockets in points 
with sharply weakened pressure, for instance 
along the line of deep breakings. The local 
decrease of the melting temperature caused 
by the weakening of pressure facilitates the 
formation of magmatic pockets (UrFEN, 1959). 

Many geophysical facts can be interpreted 
in terms of such belt that will be called the 
‘“‘ differentiation belt.’’ For instance it is 
known that the deep earthquakes are not ob- 
served at depths greater than 600 to 700 km. 
Below this depth a zone of relative stability 
begins. The lower boundary of our ** dif- 
ferentiation belt’’ is just at this depth. The 
magmatic foci of the volcanoes can be con- 
sidered as the remains of melting that have 
lasted up till now as separate cases. Actually 
according to the seismic determinations of 
G.S. Gorsukov (1956), the hearth of Klyuchev- 
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sky volcano lies under the earth’s curst in the 
upper mantle of the earth, which coincides 
with the location of the upper boundary of 
the ‘‘ differentiation belt.’’ At last at depths 
of 600 to 900km the so-called C-layer with 
high seismic velocities and a sharp jump of 
electrical conductivity is known to exit and at 
the depth of about 150km a layer with low 
velocity of seismic waves, the so-called low- 
velocity layer, is observed. These facts may 
also be explained as the traces of once existed 
molten layer that underwent the differentia- 
tion which resulted in the lifting of more acid 
rocks to the top and the lowering of the 
material more rich in iron to its base. The 
re-distribution of the radioactive elements 
caused an increase of temperature gradient 
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in the so-called low velocity layer and a de- 
crease of it at the base of the ‘‘ differentia- 
tion belt’’ which coincides with the location 
of the C-layer. Calculations show that the 
gradients are just such which are necessary 
for the observed change of the velocities of 
seismic waves (Lupimova, 1959). Thus we can 
explain the physical nature of the low velocity 
and C-layers from a new point of view. 

In accordance with the discussed theory the 
volcanic and intrusive activities were the main 
processes due to which easily fused and vola- 
tile fractions were brought to the earth’s 
surface. 

Our considerations of the thermal evolution 
of the earth leads to a conclusion that the 
volcanic and tectonic activities as well as the 


Temperature distribution with account of the variation of thermal conductivity. 
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Fig. 1. The distribution of temperature 
moment of the origin of the earth. 
process of the formation of the earth’s crust 
in general must be due to deep activity ex- 
tending to the upper layers of the mantle, 
i.e. the ‘‘differentiation belt’? mentioned 
above. 


with depth for different time, since the 


An independent consideration of the chemi- 
cal evolution of the earth recently made by 
Academician Vinocrapov (1959) draws to the 
Same conclusion. VinoGRapov conducted an 
interesting comparison of the content of H,0O, 
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CO.(C), F, Cl and other gases in the atmo- 
sphere and ocean with their content in the 
material of the earth’s mantle of dunite of 
meteorite structure. He showed that only the 
matter of the mantle could provide for the 
appearance of all the amount of volatile sub- 
stances to the surface of the earth. He gave 
many convincing facts that the source of these 
steams and gases must be a persistent vol- 
canic or intrusive activity. Vinocrapov cal- 
culated that layers to the depth of 800 to 
900km must take part in the formation of 
volatile substance. The experiments of the 
zone of melting of the chondrite ‘‘ Saratov ’’ 
showed that the result of the melting of the 
meteorite matter would be the formation of 
basaltic lava on one hand and of dunite on 
the other. It can be expected that as a result 
of such melting the maximum quantity of the 
molten light fraction would be about 6 per 
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Fig. 2. The distribution of temperature in the 
upper layers of the earth for the continental 
type of the crust for the case of quick dif- 
ferentiation of the matter ended 3-109 years ago 
(solid curves) and for the case of continuous 
differentiation (dotted curves) for the moment 
of time ({=3-0; 4-5; 5-10° years from the be- 
ginning of the existence of the earth. 


cent. This well accounts for the fact that the 
earth’s crust is rather thin and it is only 1 per 
cent of the total thickness of the mantle and 
that the matter of the mantle must transform 
into the matter of dunite type in the process 
of the evolution of the earth. 

Since the moment of the origin of mag- 
matism the volcanic activity and the thermal 
regime of the upper layers of the earth must 
have been closely connected with each other. 
The theory of geotectonic history given by 
Prof. V. V. Brtoussov (1960) emphasizes the 
significance of the carrying of heat from the 
interior of the earth in the process of the easily 
melting fraction to the surface. 

Earlier we considered the influence of the 
transfer of heat sources into the earth’s crust 
on temperature assuming that this transfer 
occurred 3-10° years ago (the age of the earth’s 
crust) practically instantaneously. This case 
is suitable for the description of the regime 
of ancient crustal regions whose formation 
ended long ago. These are ancient shields, 
possibly platforms. For this case the distribu- 
tion of temperature is shown by solid lines in 
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Fig. 3. Same as Fig. 2 for the oceanic type of 
the crust. 


20 H. A. LUBIMOVA 


Fig. 2 for the case of the continental crust 
and in Fig. 3 for the case of the 10 km oceanic 
crust. 
Another extreme case would be that the in- 
crease of the concentration of sources in the 
earth’s crust occurred during the last 3-10° 
years uniterruptedly by equal rate. In the case 
the generation of heat in granitic and basaltic 
layers of the earth increased with time, as is 
shown in Fig. 4. The heat generation in the 
mantle respectively decreased. This case was 
theoretically discussed by the author (LuBImova 
1959). Numerical calculations were performed 
by means of hydrointegrator together with Van 
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Fig. 4. Heat generation in the granitic and ba- 
saltic layers in the case of continuous carriage 
of the heat sources into the earth’s crust dur- 
ing the last 3-109 years. 


Zin Di in 1959. The continuous differentia- 
tion of the matter of the earth’s mantle was 
discussed also by Levin and Marva (1960). 
However they did not take into account that 
the transfer of heat to the surface part was 
carried out also by the masses of the matter 
heated inside the mantle. They considered the 
earth’s crust to be homogeneous as to its 
structure and its thermal conductivity equal 
to that of dunite and did not take into account 
the exciton component of the thermal conduc- 
tivity. All these factors being taken into ac- 
count we obtained temperature curves denoted 
by dashed lines in Figs. 2 and 3 for the case 
of continuous defferentiation. The amount of 
heat Q carried out by the matter that formed 
the earth’s crust is estimated by us on the as- 
sumption that all the matter of the crust is 


tranferred from the depths where the tem- 
perature was by 47=1000~3000°C greater 
than in the crust. Then Q=VioxC.dT where 
Vi is the volume of the crust and px and Cx 
are density and heat capacity of crustal rocks. 
For the period of time 4¢ this makes a con- 
tribution 4g=Q/s4t to the surface heat flow 
where s is the area of the earth’s. surface. 
Thus for the 30km thick crust at 4¢=3-10° 
years and 4T=1000°C or 3000°C, the heat: 
flow 4g=0.11 or 0.4x10-®cal/cm*sec should 
be added to the heat flow due only to the heat 
generation in the crust. This gives the value 
of the total heat flow (1.8~2.1).10-* cal/cm? sec 
for the present time. 

The temperature decrease 67 in the mantle 
from where the material is carried must be 


éT=(52 ar where V is the volume of the 


mantle. This gives 6T=150 to 300°C if we 
assume that the material is carried uniformly 
from all the depths up to1000km. This case 
is represented in Figs. 2 and 3 for the continen- 
tal and oceanic crust. In reality the degree 
of the differentiation was not uniform at all 
the depths. In lower layers it should have 
been less active than in the upper ones. If 
we assume that the maximum carriage of the 
T?K 


~ 
1800 
16007 
1400+ 2 
1200fr 
O00 F 
8007 
600} 
400 
a 
Fig. 5. The decrease of temperature in the upper 

layers of the mantle in the process of the car- 


riage of heat sources from the :depths up to 
200 km (curve 1) and up to 100km (curve 2), 


On Conditions of Magmatism Origin and Role of Volcanic Activity 21 


material was from the upper 100 to 200km 
the total temperature decrease for the time 
interval in question would be already a con- 
siderable value i.e. dT=1000° to 500°C re- 
spectively. In Fig. 5 curves designate those 
temperatures which would be observed in this 
case. 

We can see that in the upper layers of the 
earth a tendency of cooling exists while in the 
deep interior, below the level of 1000 km, 
continuous heating is under way. 

Geological data imply (BELoussov, 1960) that 
the formation of the earth’s crust was not 
uniform in space and time. 

At one and the same moment of time dif- 
ferent areas of the earth’s crust were at dif- 
ferent stages of the development. Therefore 
the discussed example of the gradual lifting 
of the materials to the upper layers of the 
earth that occurred uniformly during the last 
3 billion years cannot be considered typical 
and characteristic for the earth’s crust on the 
whole. I[t is true in some particular case. 

In reality the geological periods connected 
with the lifting of the material from the 
mantle corresponded to some intervals of time 
<3-10° years alternating with the states of 
relative rest. 

If the process of the material lifting is con- 
centrated for instance say at the interval of 
At=1-10® years then 4g=(4~10)-10-°, if 4T= 
1000 to 3000°. This can explain the anomal- 
ously great value of the heat flow observed on 
young sea ridges in the Atlantic ocean. If 
there were 2 or 3 such periods 4t of the materi- 
-~al carriage then 4g must have been equal to 
(2~5)-10-8 or (1~3):10-® each time. Then the 
heat flow must have amounted to its maximum 
several times during the history of the earth. 

This case seems to be the most probable and 
typical. The position of the flow maximum 
in time depends on the concrete history of a 
given geological region. In different areas of 
the earth’s crust the maximum of the heat 


flow could have been reached at different 
times. In the light of all that was said above 
it now seems more reasonable to put the pro- 
blem on to the thermal regime of local areas 
of the earth’s crust instead of the regime of 
the crust as a whole. This regime is differ- 
ent in different areas. Such consideration 
will depend completely on the geological 
history of the region, its past and present 
volcanic activity. 
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Abstract 


The siress-sirain relaiions for granite under various high confining 
pressures were observed experimentally. Mean Young’s modulus is nu- 
merically consianit independent of increasing pressure. The volume in- 
crease in ihe fraciure range, observed characierisiically at an atmospheric 
pressure, decays wiih pressure. The empirical formula of pressure-sirengih 
relationship is given by 


P*=P)*(kP_+1)". 


li seems that the phenomena above meniioned have ihe close connec- 
tion with compressibillty G.e. porosity). The pressure-sirengih relation 
was calculated with use of Grifhih’s crack (pore) theory, putting ihe re- 
asonable (or convenient) assumpiions inio calculaiion. This calculated 
relation, deduced from ihe empirical equation of compressibility with pres- 


sure, gives the same formula as above meniioned empirical one. 


1. In recent years, a large number of experimenial resulis on the 
deformation and fraciure of rocks have been reporied, but most of them 
were carried out for carbonaie rocks such as marble and limesione (Griggs, 
et al., 1951, Turner, et al., 1954, Roberison, 1955, Paterson, 1958). These 
rocks exhibit ihe rheological properiies the same as igneous rocks at an aimo- 
spheric pressure and room iemperaiure. Under high pressure, however, 
these rocks flow plastically showing ihe disiinct yielding zone, while the 
igneous rocks have scarcely ihe sign of plastic deformation up io the con- 
siderablly high pressure, at room temperaiure. The mechanism of ihe de- 
formaiion and fraciure of igneous rocks is much complicaied as is observed 


about britvile substances in general. As the igneous rocks, however, have 


‘ 


ihe close conneciion with the rheological behayior in ihe deeper pari of 
the crusi, ii must be made throughly clear how these rocks deform and 
fracture under high pressure. 

In the previous paper, the experimenial resulis of ihe deformaiion and 
fraciure for granite under an aimospheric pressure have been reporied 
(Maisushima, 1960). Here, we report ihe experimenial resulis under high 
confining pressure, then express a brief consideration upon ihe mechanism 
of ihe fracture phenomena. 

Of course, in order io obiain ihe definite knowledges on the earth’s 
interior, we musi iake inio consideration ihe physical and ihe chemical 
states of the earih’s interior such as iemperaiure, iis gradient and iime- 
fluctuation, heai-flow and -generaiion, siaies of ihe iniernal siresses and 
their dependence on time, and ihe consiiiuent subsiances. Bui we shall 
neglect all these effects and confine ourselves inio the siudies of deforma- 
tion and fracture of graniie. 

2. The apparaius used for this experimeni is the conveniional iriaxial 
testing cylinder as shown in Fig. 1. The shape and the size of specimen give 
the considerable affection upon ihe aspect of deformation and the sirenghih. 
Especially as to the crysialline aggregate consirucied with a large size of 
grains, such as granite, these effecis may be sirikingly large. Then ihe 
‘capacity of pressure vessel was 


taken as large as possible to be 
able to coniain ihe large size of 
specimen. Therefore, ihe durabi- 
Specimen 


lity against pressure of this vessel 


Electrode 


Contining 


was sacrificed inevitably. This ee 


vessel can endure up to 5,000 
aim., and axial siress can be 
produced by ihe 300 ion press. 
All the specimens were enclosed 
to prevent ihe confining liquid 
peneiraiing inio ihe specimens. 
Synihetic adhesive rubber was 10cm 


: : Fig. 1. Triaxial Testing Cylinder. 
used as covering maierial, for 


its good insulating characier and flexibility. The strain was measured 


by the sirain gauges of eleciric resistance type, same as in the previous 


A 


experiments. Accuracy of the measurement was within 50 aim. as to the 
confming pressure, 100 kg/cm? as io the axial compressional siress and 
10-° as to the sirain. 

3. The relations between ihe axial compressional siress and the sirain 
under various confining pressures are shown in Fig. 2. The sitains were 
measured in ihe axial direction and along the circular arc of the cylindrical 
specimen. The long columns of Kiiashirakawa biotite granite, 30 mm. in 
diameter and 60 mm. long were used up io 1800 aim. confining pressure, 
25 mm. in diameier 50 mm. long up to 3800 aim., and 20 mm. in diame- 
ter 50 mm. long above 3800 atm.. The slender specimen has a tendency 
io bend, and ihe stumpy one may deform inio barrel shape. In these ex- 
perimenis, it was scarecely observed that ihe specimens deform in barrel 
shape wiihin the elastic range. Near ihe rupiure point, however, ihe plas- 
iic flow was observed, ihough it was noi considerable, and a slight barrel 


shape deformation was recognized. 
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Fig. 2. Stress-Strain Curves from Axial Compression Tests for Kitashirakawa 
Granite. 

Each full line denotes the longitudinal direction and broken line the lateral 
direction. 

The siress-sirain curves in ihe longitudinal direciion have the approxi- 
maiely consiani slope, independenily of confining pressure, as shown in Fig. 
2. Thai is, it can be said thai the mean Young’s modulus of this rock is con- 
siant ai any confining pressure. In ihe fracture range under enough high 
pressure, however, as above meniioned, ihe curves have the iendency io be 
concave ioward the sirain axis and show the indisiinci yielding zone, pro- 
ducing ihe appreciable plasiic flow. 

In Fg. 3, the relations between ihe siress and preudo-Poisson’s ratio 


under various confining pressures are shown. Where perceniage stress is 
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used as abscissa, ihe perceniage siress is defined as a hundred times of the 
raiio of siress io the rupiure sirengih. The siress and Poisson’s raiio rela- 
tions under various pressures can be clearly shown by the use of perceniage 
siress. With ihe increase of ihe confming pressur, ihe unusual lowness of 
pseudo-Poisson’s ratio at the first siage of loading is gradually lost and 
the volume increasing effeci in the fraciure range, observed characierisiically 


ai ihe ordinally pressure (Bridgman, 1949), decays rapidly. 


a 
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Fig. 3. Pseudo-Poisson’s Ratio (ratio of lateral strain to longitudinal strain) vs. 


Percentage Stress under Various Confining Pressures. 


Rupiure sirengih (approximately equal io yield sirengih) increases sirik- 
ingly wiih the increase of pressure. The incremeni of sirengih, however, 
was noi so considerable ai fairly high pressure. ‘The observed values of 
sirengih ai various pressures are shown in Fig. 4(a). Blane circles express 
ihe values for the specimens in the ratio of lengih io diameier, 2:1, and 


full circle, 2.5: 1. 
The empirical formula which expresses the pressure-sirengih relation is, 


P* = P)*(kPgt+ 1)", 
where P* is ihe sirengih, & is the consiani and Py is ithe confining pres- 


sure. 
The imporiani numerical resulis are listed in the following Table, 
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Fig. 4. Strength vs. Confining Pressure. 

(a) Observed values and empirical curve (full line) of strength for Kita- 
shirakawa granite, and assumed strength (broken line) as the pores are perfect- 
ly closed. 

(b) Extended figure of (a), and strength of quartz. 


Table. Mean Young’s Modulus, Rupture Strength, and Indistinct Yielding Point. 


Confining | Mean Townes Modulus | puptre: Strength ige ees 
Ex10-" Stress Range kg/cm? | 5 = 
nee ee ee ee oe ieee 
Lpieea 0- 1250 1380 ° ~— 
550 | 6.9 ~ 3960 4860 4200 
10200 | 6.25 — 5940 6660 6000 
1500 6.7 ~ 7800 9100 8300 
LGOO ue ONES ~ 8530 9160 - 
2600 | 6.5 — 8580 10400 8600 
3250 | 6.25 ~11600 12110 | = 
3800 6.2 ~12500 13520 12500 
4400 6.0 10100 12150 10900 

| 


On the other hand, ihe sirengih of quariz measured by Bridgman 
(1952) is shown by mark X in Fig. 40h). The sirengih of this subsiance 
is relatively little affected by ihe pressure and only slightly increases as 
pressure does. The pressure-sirengih relation of granite at considerablly 
high pressure seems io become the same as that of quariz by the reason afier 
Grifhth’s theory (1924). Therefore, ithe assumed P*—Py curve can be ob- 
iained by exirapolaiing the above relation ioward ihe lower pressure. This 
curve is shown by the broken line in Fig. 4(a). 

We regard the difference between ihe assumed curve and ihe observed 
value as “the strength lowering’’ P*’, considering ihat the sirengih is 
lowered by a ceriain cause characierisiic of granite. 

Now, we compare ihe change of “‘the sirengih lowering” and the 
volume increase in ihe fracture siage with ihe decrease of compressibility 
(Adams, 1951). As shown in Fig. 5, there is a close conneciion among 
these phenomena. As ihe compressibility decrement with the increase of 
pressure may be caused by the closing of pores, it is reasonably consider- 
ed that ithe sirengih may be sirongly affected by the existence of pores. 
In Fig. 5, the volume increment S in ihe fracture range show ihe area of the 


part where the pseudo-Poisson’s ratio goes over 0.5 in Fig. 3 as the index. 
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Fig. 5. Comparison of ‘Strength Lowering’’, Volume Change in the Fracture 


Range and Compressibility, for Granite under Various Pressures. 


4. We will iniroduce one represeniaiion with respect to the relaiion- 


ship between the sirengih and ihe prosivy of granite under various confin- 
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ing pressures, it will be considered thai these two have a close conneciion 
with each other. 

We assume that ihe pressure is able io affect io ihe porosity alone, 
bui not to iis strength. Then the rupiure can be direcily produced only by 
the differential stress acting on ihe porous specimen. Of course, the pres- 
sure exerts the indirect effect upon ihe sirengih, changing the shape of 
the pores. 

Compressibiliiy « is given as follows, 


l Vin) CV gee, ee 
iE Va Ppae V dP htt kms Vi>Vsz, ed 


where V, Vi, Vz, the total volume 
of specimen and ihe volume of sub- 
sianiial and empty pari respeciively, 
x, denoies ihe compressibility of ma- 
terial and «por the effeciive value wiih 


respect io the pores. 


ss Put N ihe number of the pores 
2 : coniained in this specimen, and 1, 
3 Voytrr vw the volumes of ihe respec- 
2 ee oe 
ey Rockport granite tive pores, then 
(oe) 
O A P 
Quincey grante 1 “a N 
Kpr= —-=> SF (33 v:). seaw seesis @) 
V @Pa\ex J 


enclosed Ade 3 
The empirical formula of fpr is 


unenclosed : 
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Fig. 6. Compressibility of Enclosed k= “3 oe ee (3) 


and Unenclosed Specimens of Granite 


under Various Pressures (Zisman, 


from the experiment on the measure- 
1933). 


menis of the compressibility by Zisman 
(1933), considering the difference between ihe compressibility of ihe cover- 
ed specimen and ihe uncovered one. 


Now, we assume thai ihe coniraciion of ihe respeciive po-es wiih pres- 
sure has the similar tendency, then ihe effeciive compressibility of 7-ih 
pore in this specimen is given by 
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‘ ee la Bia at hee 
aes CmiPat+ni)? n®(kRPy+1)?’ | 
eA (4) 
spe pk eel | 
Ke ni ¢ in Be 
This formula (4) corresponds to the 7-ih term of formula (2), then 
pe ELT Ee) mea oes 
Ware SOs (5) 
Integrating and put the condition Paco, vs—-0 inio account, 
as ml 
Oe a EPA ED) wr Heitlatofataiciess clelecinicis's (6) 


On the other hand, we assume the shape of pore as ellipsoid rotated 


around the minor axis, then the volume is 
Anc?b; 
a 
where ¢; and 4; are the major and minor radius respeciively. 
Equaiing formula (6) and (7), and assumeing that c: does not vary with 
pressure, then the relation of the lengih of b; with pressure is given by 


pee eel oe 
Arc, myn; (kRPa+1) 
Nexi, we iry to obiain the critical siress f;, according io Grifhih’s 
theory (1921), at which the elliptic pore or crack begin to elongate by the 
applied compressional siress P along major axis, as is shown in Fig. 7. 
We iake elliptic co-ordinates a, B and 


put a=ar ai the boundary of the 7z-th 


crack. The major and minor radius weet 7 Sep 
is given by c; and 0;=a0iC;, puiiing ne 27 
: c 
c; the half length of focal line, as 2 
doz is very small. Fig. 7. Elliptic Pore Exposed under 
Then the increase of sirain en- Axial Compressional Stress P along 


ergy W; by such crack is given by eo . 


mire? 


ya 7) ie bn ee ere (9) 


W, = hot 


where E, o are the Young’s modulus and Poisson’s ratio of this material 
respectively. On the other hand, the surface energy of this crack is 


Use AceT ptt nen aneatenne (10) 
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where T is the surface energy per unit area. From the condition thai the 


erack may extend 
d =O, a. eee ll 
le, (W,—U;) =0, (11) 


the sirengih of ihis crack can be expressed as 


where @ is the numerical consianit (Ono, 1949). 

Subsiiiuiing (8) inio (12), the relation of confining pressure and ihe 
sirengih of ihe crack orienied io the direciion of compressional siress is 
given by 


1/2 
pit = az ae ann) (khPa+1)*# 


= Pit(RP bl ya, | eee (13) 
where Pio* is the sirengih of this crack at Pa=0. 
The value of & obtained from the above siated empirical equation is 


k=2.49. 10° X (kg/cm?) 7. 


Though ihe value of k for Kiiashirakawa granite from ihe compressibility 
measuremenis has not been given, the values of k from Zisman’s daia are 
k=0.22 10-?X (kg/cm?) 

for Quincy granite, and 
k=0.30 10x (kg/em?)7 
for Rockport graniie. 

5. The mechanism of deformation and fraciure of granite becomes 
much simpler at high confining pressure, as seen in ihe case of quariz 
and so on, than ai lower pressure, by the aciion of hydrosiaiic compression 
which decreases the porosiiy of porous media, while such rocks show a 
very complicated behavior under low pressure. This tendency may be kept 
at high temperature and high pressure. Therefore, in the deeper pari of 
the crust where rocks are confined with enough high pressure, the feaiure 
of mechanical behavior of rocks seems io have lost their complicacy such 
as seen at the earih’s surface, and may behave as a simple aggregate of 


constituent minerals. 


So far as it is concerned at room temperature, these soris of rocks are 


ifil 


not so much plastic even under high pressure. On the conirary, ihe elastic 


range is extended strikingly with the increase of pressure as ihe resulis of 


the raise of sirengih. However, if ihe pressure increases highly enough, the 


incremeni of sirengih does not become so remarkable, then the elevation of 


temperaiure may become more sirongly effective on the plasiicity than the 


increase of pressure, in the earih’s inierior. 
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ABSTRACT 


An equaiion of siaie of forsterite is deduced theoretically based upon 
ionic model. The calculaied variations of the density and incompressibility 
with pressure agree with those obiained by Bullen for ihe B-layer. More- 
over, the aciivaiion energy of forsieriite having Schoiiky defecis is evaluied 
io be 3.29 e.v., which is compatible with the experimenial value of 3.0 e.v. 
measured by Hughes. The variation of the aciivaiion energy wiih pressure 
is investigated, and it is resulted thai the exiraporaiion from ihe experimenial 
values obiained within 10,000 bars is noi reliable and the second derivaiive 
of aciivation energy concerned with pressure is not meaningless. The 
tempeature distribution thus obiained within ihe B-layer is well in accordance 
with Guienberg’s one. By comparing ihe laiiice energies of olivine and 
spinel (normal and inverse), which is evaluaied theoretically, we find that 


the polymorphic iransiiion from olivine io spinel is imposssible. 
1. Introduction 


The opinion thai ihe B-layer is composed of duniie is supporied widely 
with geophysical, geochemical and peirological evidences. But many problms 
are yet open io disccussion, for example, io investigated whether the variations 
of physical properiies of dunive with pressure agree or not with ihose within 
the B-layer, is one of ihese problems. 


We have only two observed quaniiiies on ihe physical siaies of maierials 
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composing ihe earih’s inierior: ihe disiribuiions of seismic wave velocities 
and eleciric conduciiviiies. Duniie is expected io have the physical. pro- 
periies deduced from ihose geophysical data. 

Olivine, main consiiiueni of duniie, coniains usually 90 per cent forsieriie 
(Mg2SiO.) and 10 per ceni fayalite (Fe2SiO,). Thereupon the physical pro- 
periies of forsieriie should play most weighiy pari in those of the B-layer. 
We will esiimaie ithe physical properiies of forsierite based upon solid-state 
physics, especially from microscopic point of view, since this kind of inyesii- 
gaiion for minerals has been little developed due io difficuliies arised from 
complicaied siruciure of minerals. 

Anoiher problem of olivine-model for ihe B-layer is concerned wiih ihe 
origin of ihe C-layer, which has peculiar physical properiies. They are an 
abnormal increase of seismic wave-velociiies and an abrupi anormaly of elec- 
iric conduciivity. Such an exiraordinary behavior of ihe C-layer yields an ait- 
raciion of many geophysisis: and enormous numbers of papers have been 
published on the origin of the C-layer. Ai present ihe following iwo types 
of hypothesis are curreni: a polymophic phase-iransiiion hypothesis which 
is realized by a change in physical siaies of silcaie, and a chemical transi- 
ition hypothesis which is realized by a chemical inhomogeneiiy there. We 
will discuss in this paper a possibility of polymophic transition of forsterite 


from olivine-siruciure io sninel siruciure. 
2. An equation of state of forsterite and its elastic properties 


According io Bragg’s result (1937), MgaSiO, crysial is orthorombic and 
has the ihree axes of ihe following lengih: | | 
a=4,.755A 6=1021A c=5,.985 A. 
In Fig.1 ihe siruciure of Mg2SiO, (olivine-siruciure) is shown. The cohe- 
sive properiy of MgzSiO, seems io be of boih ionic and covalent, and 
hence is considerably complicaied. As the first approximaiion, the assumption 
that Mg2SiO, is an ionic crysial containing Mg**, Si‘* and O?- may be -ac- 
cepiable, alihough the exact evaluation needs quaniam mechanical ireaimeni. 
As an expression for laitice energy of ionic crysial, Born and Mayer (1932) 
gave ihe following, 
u(v) = —A/(u/v)#*+ Bexp{ — (v/v)*!*/o}, D 


a 


— 5965 


—_——— = 10.214 ———___—_» 
@B O OXYGEN @O MAGNESIUM 


Fig. 1. Olivine structure of Mg,SiO,y. Si atom of which Sinema is 
the center of O-tetrahedral, is not shown in the figure. 


where v shows volume per molecule, and A,B and op are consianis but 
depend upon ihe crysial siracture and ihe constituent ions, and the suffix 0 
refers to zero pressure. We should prefer the expression (1) io another 


expression 
u(v) = —A/(v/v9)!8+ B/(v/v9)”!* 


which are also frequenily used, since (1) has been shown by Pauling 
being able to be deduced from quaniam theory, and our problem concerned 
with high pressure demandes such a formula thai is applicable io wide 
range of the variation of v/v. The first term in (1) shows elcirosiaiic poien- 
tial energy, and A is called Madelung constant, as known for such simple 
crystals as NaCl, CsCl, rutile eic., but not for olivine structure. In the 
first place, we musi evaluate the Madelung consiani for olivine siruciure. 
This evaluation is carried oui through summing up & €;/7;;(@#J), where 
é; is the charge of the tih-ion and 71; the distance between the ith-and jih ions. 
Among several practical methods for this summing up, the method deviced 
recently by Bertaut (1952) is most suitable io be applied for such a com- 
plicated erysial as olivi ine. This method is as follows: a spherically symmeiric 
charge-distribution a(r) having the center ai ihe center of each ion is re- 
placed for the actual ionic charge, where risa vector showing the position 
from the center of the ion and, for convenience, we impose on o(r) the 
restriction that each o(7) may contact, but must not overlap each other, 
that is, 


([=0 |7i>R 


where we adopt a half of the minimum distance between Si‘+ and 02 as R 


C=3/4rR', (2) 


of (2), which is the maximum value permitted without violating above res- 
tricion (2). Expanding the crystal charge density-distribution (7) 
deduced from these o(7) to a three-fold Fourier series, we obiain the po- 


tential © satisfying Poisson’s equation 


do = —4Anp(7). 3) 
The ¢ corresponds to the first term of (1), namely the electrostatic potential 
energy. In practical calculation, we adopt the cenier of Mg** as origin and 
rectangular coordinates parallel to the three axes of olivine crystal lattice 


and then summation is over all ions within the sphere having the sufficient 


radius to contain the twelfth cell (ihe unit cell is composed of 8Mg**, 45i‘* 
and 16(°~) along a-axes from the cell containing the origin. Then we 
obtain 

Ae?/(u/v9)!!8 = 155.226/b (4) 
for a molecule, where 0 is the length of b-axes. 

The second term of (1), expressing the repulsive potential energy due to 
the overlapping between each electron-cloud, has two unknown paramenters 
Band p. To deiermine these, we differentiate (1) with v once and twice, 
and then the expressions for pressure p and imcompressibility K, 

= —du(v)/dv=(1/3v)( -A/(v/v)*” 
+ BE (v/v0)1!*/o expt — (v/v0)**/o} (5) 
K= —vdp/duv=(1/9v)( —A/(v/v0)1? + BE (v/v9)!*/0 Jexp{ — (v/v9)” 
+ B{ (v/vo)?!8/o* }exp{ — (v/v0)"*/o} +P, (6) 
are obtained. For static equilibrium )=0, using the experimental value for 
the imcompressibility K of MgzSiO, at zero pressure, we can determine B 
and o. The explicit expression of (1) thus obtained is as follows : 
u(v) = 350.709185/(v/vo)1!/8-+5890.089736 expt — 4.2735(v/v0)?!*}. (7) 


From (5), (6) and (7), the variations of K, and the density d with (v/v9)'" 
are evaluated, where for (v/v)!®=1, d=3.3 gr/cm® is adopied (ihe calcula- 
tion by using aiomic weighis and lattice dimensions gives 3.216 gr/cm*), 


These resulis are shown in Table 1. Comparing these values wiih those 


obiained by Bullen (1943) for the Bayer, we can find well accordance 
between both values (see Fig. 2). In fact, both values of d agree perfecily 
and eyen the maximum discrepancy between both of K amounis to only 
8 per ceni, which may be found to correspond io difference in the seismic 
velocities for a given depth of as much as 2 or 3 per seni, as pointed out 
by Birch (1939), since unceriainiy in K/d is increased several times due 
to the fact that K/d is deduced from V.?—(3/4) Vs? (Vp and Vs are the 
velocities of compressional and shearing waves, respectively). In addition, 
considering that the calculated value denotes isothermal one and the ‘‘ob- 


served value” adiabatic, the discrepancy as above cited should be regarded 
as a slight one. 


Table 1. Variation of density d, incompressibility R and pressure p of 
Mg2SiO, versus change (v/v )1/8. 


(v/oo)i!8 | (1012 sy tes (1012 pene ees, 

1 0 | 1.22 | 3.30 | 

0.99 0.039 1,32 3.40 

0.98 0.082 | 1,44 | 3.51 | 

0.97 0.129 | 1.56 | 3.62 

0,96 0, 180 | 1.70 | 3.73 | 
— < = a a = = = | 
18 42 
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Incompressibility (10”dyens/cm’) 
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Fig. 2. Distribution of density (curve 1) and incompressibility (curve 2) 
within the B-layer. 


3. Electric conductivity of Mg.SiO, 


The eleciric conductivity of material within the Earth is deduced from 
geomagnetic daia. On the other hand, Cosier (1948) found thai the electric 


conductivity o of igneous rock is expressed by 


o =o exp(— Eu/kT) (8) 


where Ex is called activation energy, & is Bolizman’s constant and ao 
is constant. This result means that ihe eleciric conductivity of igneous 
rock increases with temperature in a similar way io that of ionic and semi- 
conduciors. Cosier, Hughes (1955) and Rikiiake (1952) estimated ihe tem- 
perature within the Earth from formula (8) in conneciion wiih the disiribu- 
ition of ihe elecirie conduciiviiies ‘‘observed’’ within the Earth. However, 
their treatments of the variation of E4 wiih pressure come inio question, 
for example, Coster’s treatment is restricted io the variation of ihe thermal- 
vibrational term of Eu, and Rikitake’s contains an unceriainiy in the 
estimation of the pressure-dependence of lattice energy and, in addition, 
disregards the variation of polarization energy with pressure (the details 
will be cited later: also see T. Wada, (1958)). At any raie, this problem 
is desired io be studied furiher. 

Hughes has shown that olivine exhibits three types of conduciions. 
Above 1,000°C the conduction becomes primarily ionic (see Runcorn and 
Tozer, (1955)). In presenti section, we may consider only ionic type of 
conduction, since ionic type prevails the other types within the B-layer, 
where the temperature is anticipated above 1,000°C. lonic conduciors are 
insulators in which certain mobile ions may move through ihe crystal lattice 
as result of defects in it. Assuming a Schoitky defect, for Me2SiO. Gn 
most of ionic crystals Schottky type is more advantageous than Frenkel 
type) (see Mott and Gurney (1948)), ihe activation energy FE is represent- 
ed with 

Fm iOl/ 2) CW aH W pW) te (9) 


where W, is the work required to remove a positive ion out of crystal, W- 
the work required to remove negative ion, Wz the lattice energy of the 
crystal per ion-pair, and U is the work required to move the ion or hole 
half way io its next metastable position. We may consider Mg** and O2- 


as the mobile ions, since Si‘t, having higher charge than Mg**, requires 


higher energy to be removed 
than Mg?*. 

(a) We can not directly 
estimate Wz, since the lattice 
energy of olivine obtained in 
the previous seciion van not 
be shared to the three kinds 
of ions. Fortunately, the elec- 
trostatic potential energy of 
spinel structure, which is very 
similar to olivine structure 


(Fig. 3 and Fig.6; see also 


section 4), has been evaluated 


QO OXYGEN @O MAGNESIUM by Verway, de Boer and van 
Fig. 3. Spinel structure of Mg,Sio4y, of which Santen (1948), whose result 

a section most similar to Fig. 1 is shown. carried out with Ewald’s method 

is separable for the three kinds of ions. Thereupon we assuem that 


buafdie= beef Lee where 
Ty : lattice ennergy of olivine structure per molecule 
Ls : elecirostatic energy of spinel structure per molecule 
Ioq : lattice energy per a-ion in olivine structure 
Isq : electrostatic energy per a-ion in spinel siruciure. 
(a= Mg?*, Sit* and O2-) 

Finally we get 


W:=35.84 e.v., (10) 


from above assumption and (7). 

(b) Ws and W- are calculated through the method deviced by Moit 
and Litileton (1937). When an ion is removed, the medium surrounding the 
vacant lattice point will find itslf in an electric field. It therefore become 
polarized, and in calculating W, or W_ the energy of this polarization 
must be taken into account. The polarization of the surrounding medium 
will set up at the vacant laitice point an electrostatic potential, which we 
denote by ¢+ or d-. Then the work necessary to remove an ion will: be, 
instead of Wz, 


Ws = Wi-(1/2)2Z.e6+ or W-= Wz-(1/2)Z_e¢-. ais) 
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Our problem is to calculate $+ and ¢-. We discuss first the problem of a 
rigid lattice, in which the ions are not allowed to move away from their 


mean positions, but in which they have polarizabilities 
amg=0.094, asi=0.0165 and a =2.4. (10-%4 cm?) 
Suppose a charge @ is placed at a definite point of the lattice, refered to 


as the point Q. Then dipoles are induced on all other ions. If 7 is the 


distance from Q, the center of the hole, we have for the polarization P 
P= (1/42){1—(//en)} (e/7*) (12) 


where ¢& is the dielectric constant for alternating fields, that is obtained 
from the polarizabilities of ions in connecton with Clausius-Mossotti’s 
formula. The dipoles on the three kinds of ions are respectively 
te={ai/ Qamg+4a0+as1) (1/42) 1 — /e0)) (ve/27*) (13) 
(¢= Mg, Si and 0) 

where v is volume per an ion-pair of Mg?* and @-. Since a dipole wz; at 
a given laitice point gives a potential y;/7* at Q, we obiain the potential 
by summation : 


B= Bmg/P+ Xtesi/P+ 2 p00/7*, (14) 


where the summations are over the lattice points of type Mg, Si and 0. 
Summations of these type have carried out with Misra’s method. We pass 
on to the case of an actual crystal in equilibrium in a siatic field, so that 
the ions are displaced into new positions. Let —cx be the restoring force, 
due to overlop forces, on any ion in a uniformly polarized medium where 
each ion is displaced by a distance +. Then dipoles g; induced at the lattice 
points of types Mg?-, Si‘* and 0% are 
pc = aif E+ (42/3) (Pot Pr) + (Zie/2e)* E+ (42/3) Pot Pr} 5) 
(¢=Mg, Si and 0) 


where E is the electric field due to absence of ion, that is, E=Ze/er?, ¢ 
is the dielectric constant for static field, Py is the polarization due to the 
induced dipoles on the ions and P, is the polarization due to the displace- 
ment of «the ‘ons. FE, ¢, Py and P; is deduced from a and c. To obtain 
directly the ¢ for Mg:SiO. is very difficult. But we can find the ¢ for 
MgO, which contains Mg*+ and 0?- and has similar elastic properties to 


Me2SiOu. (According to Szigetti (1950), we must consider the difference 
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between the ¢ duduced from the experiment of elasticity and the ¢ calcu- 


lated from the dispersion frequency.) From these valuses and (13), we. get 
(1/2)Z.e¢+ =18.63 e.v., and (1/2)Z-e$-=12.36 e.v., 


whre the ¢+ and ¢- are the potentials in the case of a deformable lattice 
corresponding to the ¢+ and ¢- in the case of a rigid lattice. 

(c) U in (9) is interpretated by Koch and Wagner as the activation 
energy of “‘structure-sensitive’’ ionic conduction, that is due to a trace of 
impurity in the crystal, acting in a way similar to the Cd ions in the silver 
salts. According to Hughes’ experiment, the conduction of olivine under 
low temperature (<600°C) is of “impurity” type, the activation energy of 
that is 1 @.v., compatible with those of common “‘structure sensitive’ con- 
ductions. Moreover, the value of oo, 10-4 ohm-! cm~}, is acceptable in 
the case of ‘‘structure sensitive’? conduction. The assumption that | e.v. 
is adopted as U may be reasonable. 


From (a), (b) and (c), the activation energy E4 becomes 


E4=3.29 e.v., 


which agrees well with the experimental value of olivine, 3.0 @.v.. 

By estimating the variations of Wz, W+, W- and U with pressure, we 
can find the variation of the activation energy E41 of MgzSiQy with pres- 
sure. (G) The variation of Wz is easily estimated, since this follows 
that of the lattice energy, given in section 2. (ii) The variation of W+ 
and W_ can not be exactly expressed, since a; and c also depend on pres- 
sure in unknown connection. We consider three cases of the variation 
of Ws or W- with pressure: (case 1) independent of pressure, (case 2) 
varies with v~1/® and (case 3) varies with v-/®, Gii) The variation of U 
also can not be exaxtly estimated, since we cannot find the repulsive 
potential in the neibourhood of the center of ion. We assume that U 
varies with v732/8, Gy) temperature also affects Eu, since ionic conduc- 
tion predominates avove 1,000°C. The influence of temperature upon E4 


is estimated by using the equation of state of Mie-Griineisen : 


p+duv)/dv=rEvp/v, (16) 
where 7 is Griineisen constat, y=1.2 for Mgz SiOs, and Eis is as follows: 
Eoiw= (9/3)@+k+3k-T-D@/T), (17) 


which shows the vibration-ierm of the lattice energy, where @ is Debye tem- 
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perature, and D(@/T) is Debye 


function. The variation of @ 


with pressure is given by @= 
@Mo(v/vo)7" (WM. =726°K for Megs- 0.010 
Si0,). The temperature effect 


TEfvib sy (1500°K) 


thus estimated is shown in Fig. 
4, where we find, at T=500, 
1,000 and 1,500°K, the thermal 
expansions of Mg:SiO; are 0.3, 


TE vib fy (1000'K) 


0.005 


0.6 and 1.0 per cent, corcespond- TE vibsy (500 KI 


du/dv (10° dynes/cm’) 


ing to decreases of pressure of 
3,000, 6,900 and 10,300 bars, res- 
pectively. 

From (i), (ii), Gii) and Gv), 10025 1005 
in the case of T=1,500°K, the (Specific volume)" 


Fig. 4. Graphical solution of the equation 
of state of Mie-Griineisen. 


variation of E4 with pressure is 


evaluated. The results are shown 


Table 2. Pressure dependence of activation energy at 1500°K. 


p | Ey (e.v.) 
ee es er) (case 1) (case 2) (case 3) 
0 3.250 oe OD 3.296 
0.010 3.290 3.290 3.290 
0.019 3.405 3.370 3,260 
0. 092 Bel 3.420 | 3,210 
0.137 | 3.620 3.480 3,150 
0.187 B20) S0h0) 3.100 


(case 1) (1/2)Zmgebug+ (1/2) Zoeh0= const 
(case 2) Y oc (y)~1/6 
(ease 3) 1 oc (y)-18 


in Table 2., from which we can find the following : d?E/dR? Ce denotes the 
distance from the center of the Earth) is not zero, as mentioned by Verho- 
ogen (1956), even in the B-layer. In fact, dE 4/dp=4.0-10-§ e.v./bar at 
p=0~10,000 bars, and dE4/dp=2.0-10° e.v./bar at about 150,000 bars. 
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Therefore, the extrapolation of dE./dp derived by any experiment under 
low pressure (/<10,000 bars) will result in overesiimating the temperature 


within the Earth’s interior, as mentioned by Hughes. 


Considering the Hughes’s value of dE 4/dp=4.3-10-° e.v./bar for peridot, 
we should take (case 1). When we adopt 10 ohm cm™ as oo and the d- 
curve given by Lahiri and Price (1939) (1) and the distribution given by 
Rikitake (1950, 1951) (2) are adopted as the distribution of the electric 
conductivities within the B-layer, the temperature distribution, is obtained, 


as given in Table 3 and shown in Fig. 4. For comparison, the temperature 


Table 3. Temperature distribution within the B-layer. 


Temperature “K 
depth (km) 
(1) (2) 
100 1,820 1,500 
200 1,920 1,540 
300 2,010 1,580 
400 2,100 1,600 
Mean of d7'/dR 1°/km 0.4°/km 


(1) calculate d value based fiper iahe and Price’s d-curve. 


(2) calculated value based upon Rikitake’s model. 


2000 


1500 


Absolute temperature (°K) 


10 
00 0 


100 200 


(7) 


Depth 


(3) 
ne 


400 500 
(km.) 


Fig. 5. Temperature distribution within the B-layer. 
(1) calculated value based upon Lahiri and Price’s d-curve. 
(2) calculated value based upon Rikitake’s model. 


(3) gutenberg’s result. 
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distribution given by Gutenberg is also shown in Fig. 4. 
4. Possibility of the polymorphic transision of Mg.SiO, 


Bernal suggested first the possibility of polymorphic transition of olivine 
at meeting when Jeffreys reported the existence of the C-layer. The 
magnesium germanate (Mg:GeQ,) has a similar chemical formula to Mg:~ 
SiO. and shows a polymorphic transision. Thus he infered that the same 
kind of polymorphic transition might take place for Mg2SiOs. Goldschmidt 
reported that Mg2GeQ, has two types of possible crystal structures, viz. 
olivine type (orthorombic) and spinel type (cubic). Many other experi- 
menial works have been carried out to confirm the above report. Romeijin 
and Ringwood found that MgzGeO, changes its structure from spinel io 
olivine at about 600°C. On the other hand, Bridgeman found by high 
pressure experimenis that the olivine shows a discontinuous jump of shear- 
ing strengih at 84 kbars. He suggested that this discontinuous Jump may 
be related to its polymorphic transition of crystal structure. 

A theoretical study of the polymorphic transition of crystal has been 
carried out on ionic crystals, especially on a simple crystale such as alkali 
halides. One way of atiack is io calculate the lattice energy based upon a 
classical treatments. Though the results thus obtained give transiiion-pres- 
sures, the values are so unreliable that only the possibilities of transitions 
are confirmed, even if higher terms of energies are considered and ihe 
dependence on pressure is estimated as exactly as possible. (see, for 
instance, Jacobs (1938), and Born and Hung (1954)). 

A stability condition of crystal at high pressure and at absoluie zero 
temperature is given by examining the Gibbs free energy (or enthalpy) 
U+PV. U,P and V denoie lattice energy, pressure and volume. When 


ihe transition iakes place at the pressure P, the stability condition gives 
UstPVs=Ust+PVe 18) 
or 
P=(Uz—U.4)/(Va— Va) (19) 
Neglecting the effect of pressure upon U-and V for the first approxima- 


tion (as above meniioned, we should be concerned with the possibility of 


transition only), we get 
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P= (Up,0— U4.) /( Vaso — Vz,0) (19') 
where the suffix 0 denotes the value at P=0. 

The repulsive energy represented by exponential type in (1) is a 
rigorous expression. However the inverse power type is frequencly used 
for many practical problem owing to its simple expression. In the fol- 
lowing we use the lattice energy as 

U= —(A/r)+ (B/r) (20) 

From the equilibrium condition (dU/dr)=0, the value 7% and Up at 

P=0 are given by 
7 = (nB/A)!®), Uy= ~(A/10)/{(n—1)/n} (21) 
From (19’) and (1), we get 


Tetrahedral 


(c) 
Fig. 6. The structure of spinel. 

The structure in idealized form. The oxygen atoms (shown by large circles) 
are approximately in cubic closed packing. (a) and (b) show two close-’ 
packed layers projected upon the cubic face. One set of silicon atoms (shown 
by small circles A) is in positions of four-coordination between a tetrahedral 
group of oxygen atoms and the other set of magnesium atoms (shown by 
small circles B) in six coordination between an octahedral group. (shown 
in (c)) We may picture as a series of layers as above shown. The bottom 
layer of (a) has diagonal chains of magnesium-octahedra,. which are linked 
laterally by silicon-tetrahedra lying alternately above and below heavily 
outlined oxygen atoms. In the next layer (b), the slant of the chains is 
reversed. Four such layers makes up the complete unit cell. 
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P=(A/sr) Fa (22) 
Fa={(—1)/n}{1—(Ax/Aa oY 
(Ba/Bz) }/{(1— (Sa/s.4) (Ba/Ba)*! (A 4—Ap)9!™- } 
where Va=sara?, Va=Sarz'. 

The spinel XY20, has generally the spinel structure which is shown 
in Fig. 3 or Fig. 6 for the normal spinel structure. We notice that ihe 
Fig. 3 is quite similar to the olivine structure of XY20, shown in Fig. 1. 
In the normal spinel X-ions occupy the center of ihe O-tetrahedra (half of 
Y-ion in the inverse case). According to Verwey et al. the elecirostaiic 
potential energy of the spinel expressed by 

Usp= — 1/2) (€?/@) (b Meetrat+29Moct+8My) = —M(e?/a) (24) 
where 
a : length of the cube edge of the cubic lattice 
Metra: coniribution from tetrahedral interstices 


Moecta : contribution from octahedral interstices 


M, _ : coniribuiion from O lattice points 
p : ionic charge number of the tetrahedral interstices 
qd : ionic charge number of the octahedral interstices 


For the normal and inverse spinels p=4, 7=2, and 2, 1/2) (4+2) =3, 
respeciively. Table 4 is the numerical iable for M. According to Bron 
and Mayer’s expression, a repulsive potential Ze between two kinds fo ions, 


a and B, 


Table 4. Electrostatic energy of normal and inverse spinel for several values 
of oxygen parameter (unit : (¢2/a) 10-4 erg) (after Verwey et al.) 


pee a Me. Mone M, i M : 
De eee) tree | Aaier | 38,6 
cS a ae eee 
ja [32 | 28 | Bo] Re | me 
cw} ap | ne) ee | Re | Be 
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R=Cha,pbexp{ Ya +1p —%axp)/0} (25) 
where C and b are the coordination number and a constant, respeciively, 
and Aa,g is expressed by » 

hap =1+Za/Nat+Zp/Ne (26) 
where Z is the charge of ion, N is the number of electron in the outer- 
shell of ion, and the suffix a or 8 denotes a- or f-ion, 7a, 7g and YasB in 
(25) are basic radii of a- and $- ions and the nearest distance between both 
ions, respectively. _ Refering to (25) and (26) we obiain the repulsive 
potential R per molecule as follows : 


for olivine-siructure 


Roe = 4dsiobexpl (7s¢+ 1% — 0.17770) /p} 
+ 12g, obexp{ (7an+ 7 — 0.20445) /p}, (27) 


for normal spinel-siructure, 


ne= 4dsiobexp{ (rsit+r—V 3 .a/8)/0} 
+12amgobexp{ (7ug+%—a/4)/0}, (28) 
and for inverse spinel-airucture, 
IS 4am obexp{ Crug +o a V 3 a/8)/o} F Amo .obexpt (mg +%o— a/4)/o} 
+ 6Asi,obexp{ (7si+%—a/4)/o}, (29) 
where the constants p and 6 are 0.345A and 10-1? erg, respectively. 

To compare the energies of the three types of structures, viz. olivine, 
normal spinel and inverse spinel, it is convenient to transform the unit of 
length from @ (cubic edge length of spinel) to 6 (b-axes length of olivine), 
By using this transformation, we set up approximately 

exp(—0.17776/o)/exp( — 0.20446/o) 
=exp(— 0.1770b/0)/exp( — 0.204346/p)=G (30) 

Then, referring to (22) and (30), B4/Bz is expressed approximaiely 
by Roe, Rns and Ris. 

We musi consider the effect of oxygen parameter, as previously men- 
tioned. It is well known that the actual spinel has w greater than 0.375. 
In the present section, the effect of oxygen parameter greater than = 0.375 
is discuussed. Putting w=0.375+0, then we replace VY 3a/8 and a/4 in 
(26) and (27) with 1/ 3 a{(1/8)—4} and a{(1/4) - 6}, respectively. 

In conclusion we have Fn and P in (23) as shown in Table 5, where 


35nS7. The condiiion of transition from A to B is given by Fn>0, and 
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Table a Signs of F,, and Vor/Vsp for several values of » and uw. The 
sufficient condition for a polymorphic transition is F,,>0 and Vor/Vep>1. 


u=0,375 | u=3,80 
n == = 
normal spinel inverse spinel | normal spinel __ inverse spinel 
3 negative | negative negative negative 
4 negative positive negative negative 
| <1 
5 negative positive negative negative 
on ee 
6 negative | po pak negative negative 
| ei 
i negative positive negative negative 
<1 


Vu/Vr=1. In iable5 Vu/Vz is shown only for F,>0. The results thus 
obtained mean thai ihe polymorphic iransiiion from olivine io spinel (nor- 


mal or inverse) is impossible. 
5. Conclusion 


The variations of ihe elasiic and elecirical properiies of MgzSiO. with 
pressure agree well with those deduced from seismic daia and geomagneiic 
daia within the B-layer. It seems true that the B-layer is composed of 
duniie. 

On the other hand, ihe possibility of polymorphic transition from 
olivine to spinel is denied at any raie. However io confirm ihe possibility 
of phase-iransiiion theoretically, generally speaking, is a difficult problem, 
since wheiher ihe iransisiiion is possible or noi, depends upon the difference 
between the laiiice energies of both phases. Even if the lattice energies of 
both phases is considerably well evaluated, the difference may be not reliable, 
though the unknown errors in both lattice energies are somewhat cancelled. 


Thereupon we should say thai this problem is yei open io discussion. 
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